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I. INTRODUCTION 


During the past three decades, enormous amounts of resources have been expended in the 

• design and development of Liquid Oxygen/Hydrocarbon and Hydrogen (LOX/HC and LOX/H 2 ) 
rocket engines. A significant portion of these resources have been used to develop and demon- 

• strate the performance and combustion stability for each new engine. During these efforts, many 
analytical and empirical models have been developed that characterize design parameters and 
combustion processes that influence performance and stability. Many of these models are suit- 
able as design tools, but they have not been assembled into an industry-wide usable analytical 
design methodology. 

An objective of this program was to assemble existing performance and combustion sta- 
bility models into a usable methodology capable of producing high performing and stable 
LOX/Hydrocarbon and LOX/Hydrogen propellant booster engines. In the first phase of this pro- 
gram, we focused on developing an analytical design methodology to perform the design trade- 
offs leading to an optimized high-performing and stable combustion device. This design 
methodology, the Rocket Combustor Interactive Design (ROCCID), contains previously devel- 
oped and available analysis models for characterizing the effects of all critical design and oper- 
ating parameters on the performance and combustion stability of liquid propellant combustors. 

This methodology was released for industry use in May 1991 (Ref. 1) after undergoing a 
Beta test by potential users, including propulsion contractors, government facilities and 
interested Universities. 

Concurrently, a validation plan for demonstration of the ROCCID capabilities and identifi- 
cation of its shortfalls and limitations was developed. The methodology was applied to develop 
a combustor design and a detailed test plan was formulated to verify the performance and com- 
bustion stability of the combustor compared to that predicted by the ROCCID Methodology. 

The validation hardware consisted of a thrust chamber for operation with LOX/RP- 1 pro- 

• pellants and containing an injector assembly with 105 OFO triplet elements and a combustion 
chamber with a silica phenolic liner. The thrust chamber was designed for operation up to 

' chamber pressures of 1,750 psia and over a wide range of mixture ratio (e.g., 1.0 to 10.0). A 

copper insert was provided to permit operation with and without an acoustic resonator. The basic 
hardware design is fashioned after a test proven thrust chamber developed and hot fired during a 
successful program conducted by Aerojet for the Air Force Astronautics Laboratory (Contract 
F0461 1-85-C-0100). 
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I. Introduction, (cont) 


Hot fire tests were performed with this thrust chamber resulting in 27 valid tests for direct 
comparison with the ROCCID methodology predictions. Test results and comparison with the 
ROCCID based predictions are included in this report. 

This program was conducted by the Aerojet Propulsion Division under Contract NAS 3- 
25556 from the NASA Lewis Research Center. This program was conducted over the period 
from December 1988 through November 1991. Mark Klem was the NASA program manager 
during this entire effort. 
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II. EXECUTIVE SUMMARY 


Up to this time, there has been no industry standard methodology to aid in the design and 
analysis of combustion stable, high performance liquid propellant rocket engine combustion 
chambers. While NASA design guidelines for injectors, chambers, and stability devices (NASA 
SP-8089, 8087, and 8113, respectively) were available, designers used their own methodology 
for developing design features and evaluating effects of design or operating parameter changes 
on combustion stability and efficiency. This was usually a time consuming process, requiring 
extensive designer and analysis specialist interaction, and generally dependent upon expensive 
experimental evaluations with many design iterations. 

This problem was addressed during this effort through the development of a user friendly 
and interactive computer based analytical design methodology for creation of liquid propellant 
combustor designs that produce both high performance and stable operation. This ROCket 
Combustor Interactive Design (ROCCID) methodology consists of an assembly of existing 
industry-wide performance and combustion models linked with an interactive computer logic 
that allows the user to design a combustor to meet desired performance levels and combustion 
stability margins. 

The development of ROCCID provides a design and analysis tool to reduce future engine 
development costs. Propulsion contractors and university research centers benefit through access 
to the best available analysis codes and improved design efficiency from the use of properly 
linked and adequately documented analysis tools. Efficiency will evolve from a standardization 
of the combustor design approach. Government agencies also benefit by having a valid tool for 
evaluation of proposed designs, realistic assessment of performance goals, and effective system 
optimization early in a development program. ROCCID is already in use in support of the 
STME engine for the NLS where the consortium formed for that development program has 
recommended its use. 

Methodology Overview 

The structure of ROCCID is illustrated in Figure 1. ROCCID contains three main 
components which are: 

1. An interactive front end (IFE) that provides guidance to the user for input setup, input 
and output control and the generation and maintenance of library files for replay and 
restart, propellant properties and combustion gas properties. 
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Figure 1. Major Components of ROCCID Include the Interactive Front End, 
The Point Analysis Evaluation and the Point Design Evaluation 
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II. Executive Summary, (cont) 


2. A point analysis option that provides performance and combustion stability analysis 
of existing combustor designs. 

3. A point design option that creates the essential combustor design features for a high 
performance and stable rocket engine from specified requirements. 

Both the point analysis and point design options use a library of performance and com- 
bustion analysis models selected from an existing and industry-wide inventory. These analysis 
models are contained within ROCCID in a modular format. This permits the user to access spe- 
cific models for a specialized sub- analysis or to use two or more models that perform similar 
functions to define and resolve uncertainties in the particular area of the analysis. Modular con- 
struction also permits easier and less costly methodology upgrading as new analysis models are 
developed or refined. 

A steady-state combustion analysis including propellant atomization, vaporization and 
mixing supplies key input into the performance and stability analyses. Specific models that were 
reviewed and selected for use in the ROCCID steady-state combustion analysis are described in 
Table 1. 

Four models/correlations for propellant dropsize are included for standard injector ele- 
ments from showerhead, doublet and triplet impinging elements and shear and swirl concentric 
tube (coaxial) elements. Dropsizes from all applicable correlations are calculated and displayed 
for comparison. A user may select any of the calculated values for the steady state combustion 
analysis or provide another estimated value in their place. Propellant (fuel or oxidizer) 
vaporization is calculated using the generalized length correlation developed by Priem and 
Heidmann. Propellant mixing is based on the use of a unielement mixing efficiency value 
determined from cold flow measurements and adjusted for interelement mixing effects. 

Currently, this value is supplied by the user, but guidelines for its selection are included. 

The combustion stability analysis is made with a large array of models used to calculate 
the chamber admittance and the burning response magnitudes. These models provide the 
capability to estimate combustion stability margin for all common types of combustor 
instabilities including chugging and chamber acoustic coupled (high frequency) modes. 

ROCCID displays the stability analysis results in terms of the calculated growth coefficient (A.) 
for the particular acoustic mode of concern. This growth coefficient represents the amount 
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TABLE 1. EXISTING STEADY-STATE COMBUSTION MODELS WERE EVALUATED FOR USE IN ROCCID 


I. Atomization 


Model/ 

Correlation 

Developed By 
Developed For 

Approach 

Applicable Injection 
Element Type{s) 

Advantages/Disadvantages 

Used in 
ROCCID 

Priem TR-67 

NASA-LeRC 

In-House 

Derived Empirically From 
LOX/Heptane Tests 

Showerhead, Doublet 
Triplet 

Propellant Properties Effects 
Included, Historical Data Base. 
Limited Off-Design Capability 

Yes, as Option 

Aerojet 

Impingement 

Model 

Aerojet 

NASA-LeRC 

Potential Flow/Boundary 
Layer Breakup 
Calculation 

Doublet Triplet 

Mechanistic, Simple Off-Design 
Capability Total Time Lags 
Calculated 

Yes, as Option 

PLC 

Aerojet 
NASA MSFC 

Potential Flow/Boundary 
Layer Breakup for 
Gas/Liquid Concentric 
Tube Elements 

Shear Coaxial; Swirl Coaxial 

Mechanistic, but Little Historical 
Basis 

Yes, as Option 

CICM 

Rocketdyne 
NASA MSFC 

Uses Empirical 
Coefficients for LOX 
Breakup and Intra- 
Element Mixing 

Shear Coaxial 

JANNAF Reference but Difficult to 
Incorporate and Limited Historical 
Base 

No 

SDER 

Rocketdyne 
Phillips Laboratory 

Uses Empirical Dropsize 
Correlations From Hot 
Wax Experiments 

Doublet Triplets 

JANNAF Reference but Contains 
Known Correlation Problems 

No 

Dropmix 

WJSA 

Phillips Laboratory 

Empirical Dropsize 
Correlations 

Doublet, Triplet; Shear 
Coaxial 

Improved Correlations Over SDER 

Yes, as Option 


II. Vaporization 


Model 

Developed By 
Developed For 

Approach 

Applicable Propellants & 
Operating Conditions 

Advantages/Disadvantages 

Used in 
ROCCID 

SDER 

Rocketdyne 
Phillips Laboratory 

Simultaneous Solution of 
2-Phase Equations for 
Gas Flow and Drop 
Acceleration and 
Vaporization 

All Hydrocarbons with 
Definable Properties. Sub- 
Critical and Super-Critical 
Pressures 

Contains Fatal Calculation Problems 

No 

PCDER 

WJSA 

Phillips Laboratory 

Same as SDER 

Same As Above 

Improved Version of SDER-Runs on 
PC 

No, but Could 
Be Used 
Outside 
ROCCID 

Generalized 

Length 

NASA-LeRC 

In-House 

Simplified Correlation of 
Mechanistic Vaporization 
Model Calculations 

Derived from Subcritical 
Pressure Results but Has 
Been Applied Successfully 
to Supercritical Conditions 

Simple, Fast Calculation; Excellent 
Historical Basis 

Yes 


III. Mixing 


Model 

Developed By 
Developed For 

Approach 

Applicable Injection 
Elements 

Advantages/Disadvantages 

Used in 
ROCCID 

LISP SDER 

Rocketdyne 
Phillips Laboratory 

Uses Spray Coefficients 
Derived from Unielement 
Cold Flow Tests 

Doublet and Triplet Built-In 
Coefficients 

Spray Coefficients Derived From 
Cavitating Flow Invalid, Time- 
Consuming Set-Up 

No 

Dropmix 

WJSA 

Phillips Laboratory 

Uses Intra-Element 
Mixing Efficiency Based 
on Cold Flow Data Base 

All for which Em Correlation 
Is Available 

Simple, Fast Calculation 

Yes, as Option 

Two Stream 
Based on Em 

Aerojet 
NASA MSFC 

Uses Intra-Element or 
Inter-Element Em Based 
on Cold Flow or Hot Fire 

All for which Em Correlation 
Is Available 

Simple, Fast Calculation 

Yes, as Option 

Hersch 

Turbulent 

Mixing 

NASA-LeRC 

In-House 

Based on Turbulence 
Intensity 

All 

Simple, but Does Not Account for 
Intra-Element and Combustion 
Effects 

No 
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II. Executive Summary, (corn) 

of amplification required by the chamber to achieve the condition where the driving required to 
support the waveform exactly equals the driving response present within the system. Thus a 
value of X = 0 represents a neutral stability condition, X>0 represents an unstable condition and 
X<0 represents a stable operating condition. 

The effects of damping devices such as acoustic cavities (1/4 wave tube and helmholtz) 
and baffles are also considered by these models. A listing of the combustion stability models 
considered for ROCCID is provided in Table 2. 

The performance of the combustor is defined by its energy release efficiency. This 
accounts for combustion efficiency limitations resulting from incomplete propellant vaporization 
and/or mixing. The energy release efficiency is calculated using the JANNAF simplified per- 
formance procedure as contained in CPIA 246 and the propellant vaporization and mixing limi- 
tations from the steady-state combustion analysis. An input file for use in the TDK/BLM com- 
puter analysis is also generated that can be used for a rigorous performance analysis of the com- 
plete rocket engine as described in CPIA 246. 

The design requirements for combustor cooling are established using techniques outside of 
ROCCID. These requirements may include estimates of fuel film cooling required for chamber 
and baffle walls, dump cooling off baffle tips, and bulk temperature increases resulting from 
regenerative cooling of the nozzle chamber and resonator/baffle components. This information 
is used to calculate the propellant injection temperatures, injection orifice requirements and the 
local flow injection mixture ratios. This method of accounting for temperature limits of the 
injector/thrust chamber materials with ROCCID was selected to keep the focus on the com- 
bustion stability and performance issues while providing a useful and practical combustor design 
tool. 


Capabilities and Limitations 

ROCCID has been specifically formulated to be applicable to combustor designs for 
LOX/HC, LOX H2, and N2O4/MMH propellants. Propellant and combustion gas properties for 
LOX/RP-1, LOX/CH4, LOX/C3H8, LOX/H2, and N2O4/MMH are included. Both gas/liquid 
and liquid/liquid propellant injection are considered. Conventional impinging like doublet and 
triplet (OFO and FOF) elements, impinging unlike double elements, non-impinging showerhead 
and coaxial tube elements with or without liquid stream swirl are treated. Injection element 
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TABLE 2. COMBUSTION MODELS ARE INCLUDED FOR ALL ASPECTS OF STABILITY MODELING 


L CHAMBER RESPONSE 


Model 

Developed By 
Developed For 

Approach 

Applicable Design 
Features 

Advantages/ Disadvantages 

Used in 
ROCCID 

HIFI 

Aerojet 

Phillips Laboratory 

Linear Perturbation 
Technique With Mean and 
Fluctuating Components 
for Dependent Gas 
Dynamic Variables 

Acoustic Resonators 

Mechanistic, Burning Rate and 
Injection Coupled Extensive 
Application History 

Yes 

DIST30 

Colo State 
Phillips Laboratory 

Calculates Baffle Damping 
Using a Turbulent 
Boundary Layer Model for 
Viscous Dissipation 

Baffle Height and Blade 
Distribution Acoustic 
Resonators as 
Secondary Damping 

Distributed Combustion, 
Mechanistic, Radial Baffles Only 

Yes 

FDORC 

Colo State 
Phillips Laboratory 

Piecewise Distributed 
Combustion With Arbitrarily 
Located Resonators 

1/4 Wave Tube and 
Helmholtz Resonators 
and Liners 

Distributed Combustion, Resonator 
Location, Mechanistic 

Yes 


II. BURNING RESPONSE 


Model/ 

Correlation 

Developed By 
Developed For 

Approach 

Applicable Injector 
Types 

Advantages/ Disadvantages 

Used in 
ROCCID 

Reardon-Smith 
N/t Correlations 

Reardon-Smith 

JANNAF 

Correlation of Empirical N/x 
Using a Sensitive Time Lag 
Model 

Doublets, Triplets; 
Coaxial 

Simple Historical Data Base, 
Non-Mechanistic 

Yes, as Option 

CRP 

Aerojet 

Phillips Laboratory 

Uses Agosta-Hammer 
Non-Unear Vaporization 
Response Model 

All for Which a 
Representative 
Dropsize Exists 

Mechanistic, but Can Require Long 
Run Times 

Yes, as Option 

Empirical/Damp 
Growth Rate 
Correlation 

Aerojet 

Phillips Laboratory 

Use Observed Damp or 
Growth Rates to Infer 
Combustion Response 

All for Which Empirical 
Growth or Damp Rates 
Exist 

Requires Experimental Data Base 
but Is a Means for Anchoring 
Stability Model 

Yes, as a Means 
for Increasing N/x 
Data Base 


III. INJECTOR RESPONSE 


Model/ 

Correlation 

Developed By 
Developed For 

Approach 

Applicable Injector 
Types 

Advantages/ Disadvantages 

Used in 
ROCCID 

LFCS 

Aerojet 

Phillips Laboratory 

Wenzel & Szuch Lumped 
Parameter 

All with Definable Total 
Timelag 

Simple, Good Track Record for 
Chug Instability 

No 

INJ 

Aerojet 

NASA-LeRC 

Lumped Parameter 
Analysis with Spacially 
Varying Acoustic Wave in 
the Chamber 

All with Definable Total 
Timelag 

Computes Injector Response 
Based on Element Timelag. 
Dependent on Good Injection 
Model 

Yes 

Lewis Non- 
Linear 

NASA-LeRC 

In-House 

Modification of Feiler and 
Heldmann Feed System 
Coupled Instability Model 
to Include Manifold 
Acoustic Effects 

Concentric Tube 
Elements 

Include Flow Response Due to 
Manifold Acoustics if Important 

Yes 


log 891.255A 
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II. Executive Summary, (cont) 


zones for the injector core, barrier and fuel film cooling can be included. For the point design 
option, design parameters such as contraction ratio (CR) and chamber length (L') are internally 
defined to provide the best tradeoff between performance and combustion stability. A simple 
tradeoff between nozzle length and chamber length is also included to optimize engine delivered 
specific impulse for an envelope limited system. Acoustic damping devices are also 
recommended and their design features specified to provide the required combustion stability 
margin. Design tradeoffs for a throttling engine are also performed. 

The point design option features an optimization of the injector element design through 
user interactive operation. Guidelines are provided to aid in injector element selection for a par- 
ticular application. The quantity of elements and the injector orifice size are calculated through a 
series of trade studies to satisfy: (1) performance and high frequency combustion stability 
requirements, (2) chug stability and pressure drop constraints including throttling requirements 
and (3) basic injection area and combustion chamber size limitations. 

This ROCCID methodology has been prepared with certain limitations in order to simplify 
its construction and guarantee its timely and affordable development. For example, no 
supersonic nozzle effects are included. Nozzle design and engine specific impulse are deter- 
mined outside of ROCCID. Precombusted fuel (stage combustion cycle) is not presently con- 
sidered. Mass addition from ablation, igniters, gas generator dump or transpiration cooling is not 
modeled. A flat faced injector with one element type in each zone (core, barrier, etc.) is 
assumed. As previously noted, cooling requirements are defined outside ROCCID, but their 
effects on performance and stability are considered. 

The computer code for ROCCID is operational on VAX 8600 series computers at both 
Aerojet Propulsion Division in Sacramento and NASA Lewis Research Center in Cleveland. 

This program was also operational on a SUN computer at SEA, Inc., in Carson City. Interactive 
graphics for a Tektronics 41 XX terminal are also provided. The code has been constructed 
without machine dependent instructions (except for graphics), but operation on other systems 
remains to be verified. 

ROCCID Demonstration and Verification 


The capabilities of ROCCID were demonstrated by using this newly developed 
methodology to create several combustor configurations some of which were subsequently 
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II. Executive Summary, (cont) 


designed, fabricated and test fired. These tests results were compared to the performance and 
combustion stability predictions made by ROCCED during the design formulation. This 
comparison provided a general assessment of ROCCID capabilities and a means to identify the 
strength and weaknesses of the methodology. 

The propellant combination of LOX/RP- 1 was selected for study using ROCCID and for 
eventual hot fire testing to validate ROCCID capabilities. This propellant was selected because 
historically, it has proven to be a difficult propellant to provide both stable combustion and high 
performance. 

Normal design propellant injection temperatures of 174.4°R for liquid oxygen and ambient 
condition (515°R) for RP-1 were selected. These conditions represent normal propellant storage 
temperatures within the Aerojet test facility and therefore are the most cost effective test 
conditions. 

A nominal design mixture ratio of 2.8 was selected since this approximates an optimum 
performance design from a thrust chamber performance viewpoint. Note that this mixture ratio 
value was selected for design purposes only. The effects of mixture ratio variation on a fixed 
design were evaluated using ROCCID and mixture ratio was a key operating parameter for 
validation testing. 

Thrust chamber size is a key design parameter, particularly from a combustion stability 
standpoint. The chamber diameter has a direct effect on the resonant frequency of acoustic 
modes within the chamber. Large chambers appropriate for booster engine application (Dc = 

17.5 to 44 in.), have lower resonant frequencies that are more likely to couple with typical com- 
bustion responses. On the other hand, subscale chamber diameters of 5.5 and 7.68 in. have been 
successfully used in the past for combustion stability investigations (see Section IV,C.). For this 
study, a nominal baseline chamber diameter of 7.68 in. was selected since it is a proven 
"subscale" design and residual hardware was available. 

A contraction ratio of 2:1 that corresponds to a throat diameter of 5.43 in. and nominal design 
chamber pressure of 1250 psi were selected to provide a large operating chamber pressure 
variation within the existing hardware and test facility capabilities. The operational envelope of 
the validation hardware within the Aerojet E-4 test facility is shown in Figure 2. A variation in 
mixture ratio from approximately 1 to 10 was possible at the nominal chamber pressure of 1250 
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Figure 2. Facility Operating Limitations 
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II. Executive Summary, (cont) 


psi. A maximum chamber pressure of approximately 1800 psi was achievable at the nominal 
mixture ratio. The mixture ratio excursion range diminishes as chamber pressure is increased 
because of the facility run tank pressure limits. 

Impinging injector designs were selected for this evaluation since they are most appro- 
priate for liquid/liquid injection. Specifically, the like doublet and OFO triplet elements, both of 
which are included in ROCCID, were evaluated. The OFO triplet element was selected over the 
FOF since the propellant density and mixture ratio will result in equal orifice diameters with an 
OFO element but not with the FOF triplet. From experience, equal orifice diameters are desir- 
able with a triplet element to provide optimum mixing and atomization. 

Three injector designs, fine triplet, coarse triplet, and like-on-like doublet, were selected to 
be evaluated by ROCCID for detail design and testing. Initially, the point design option of 
ROCCID was used to size the injectors for the nominal operating conditions. Thereafter, the 
analysis portion of ROCCID was used to perform trade studies on the initial designs. The 
experience gained from past experimental results were also applied to the analysis. Also from 
past experience, certain model combinations in ROCCID were used because they produced good 
results. HIFI was used to analyze the chamber response. Combination of Smith-Readon and 
Aerojet N/x correlations were used to analyze the burning response. INJ was used to analyze the 
injector response. 

The injector orifice sizes (diameter) evaluated were 0.090 in. for the fine triplet, 0.159 in. 
for the coarse triplet, and 0.100 in. LOX and 0.065 in. fuel for the like-on-like doublet. The pre- 
dicted performance efficiencies for all three injectors were greater than 97 percent so perfor- 
mance did not become a critical concern in the selection process. 

Criteria were established to select the best injector design for use in the validation of 
ROCCID. Because of limited resources, only one thrust chamber design could be selected for 
design and testing. The criteria were established on the basis of providing the best test of the 
capabilities of ROCCID. The criteria are provided in Table 3. 

The chug stability limits of the three injectors were defined using ROCCID. The fine 
triplet and the doublet have low chug thresholds at about 300 psi. The coarse triplet has a more 
undesirable chug threshold at about 600 psi. Considered that the nominal operating pressure is 
1250 psi, this higher chug pressure limit for the coarse triplet does not allow for much leeway in 
throttling the engine to find stable and unstable high frequency stability regions. 

The high frequency stability characteristics of the three injectors were also predicted using 
ROCCID. The coarse triplet was found to be the most stable injector design without acoustic 
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TABLE 3 


INJECTORS COMPARED AGAINST RATING CRITERIA 


Rating Criteria 

Injector Configuration 

Fine O-F-O 
Triplet 

Coarse O-F-O 
Triplet 

Like- on- like; 
One Doublet 

Largest negative and positive 
growth coefficients for least 
damped mode without acoustic 
damping 

Host unstable 

Host stable 

Unstable 

Best sensitivity of combustion 
stability to P c and 0/F 

Significant variation 
with 2T and IT; stable 
operation with 0/F and 
P c 

Huch larger 
variation in P c 
and 0/F required 
to drive unstable 

Good sensitivity with 
0/F for IT mode 

Greatest change from existing 
data base 

N„ 39 to 105 
D,,, 0.125 to 0.090 

N„ 39 to 48 
D„, 0.125 to 0.159 

N„ 105 to 123 
D 0 , 0.059 to 0.065 

Capable of stabilization at 
nominal operating conditions 
using an acoustic cavity 

Three types of 
cavities evaluated 

Stable without an 
acoustic cavity 

IT monotune 
cavity judged to 
be sufficient 

Best confidence in modeling 
drop size and mixing 
efficiency 

r m reasonable; 
Em good 

r m reasonable; 
Em good 

r m reasonable; 
E m reasonable 

Ease of incorporating design 
features into existing hard- 
ware 

New faceplate only 

New faceplate and 
injector core 

New faceplate and 
injector core 

Preferred rating 

6 of 6 

3 of 6 

2 of 6 



Indicates preferred 
characteristic 
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damping devices. From a flight hardware designer's point of view, this would be an excellent 
design. However, since our validation criteria require both stable and unstable operation for ver- 
ification, the coarse triplet injector design was not sufficient for ROCCID validation. If 
resources were available, this design would be a good choice for a secondary, very stable, 
alternate test series to demonstrate ROCCID's capability to stabilize the thrust chamber through 
injector design changes (i.e., combustion response). 

The doublet was predicted to be stable except for the first tangential mode. The doublet 
looked very promising because it was predicted to have a region of operation where it could be 
driven stable and unstable by changing the mixture ratio. The doublet was also predicted to be 
stable using acoustic damping devices, which would allow validation of the models for damping 
devices. The major drawback of the doublet was that it is very similar to previously tested 
designs. Since this was a validation of the predictive capability of the models in ROCCID, using 
an injector so close to previously tested hardware made this injector a secondary choice. 

The fine triplet design was predicted to be unstable in several different modes with varying 
operating conditions. Figure 3 shows ROCCID predictions of stability changes over the oper- 
ating map. The first longitudinal (1L) mode is stable until lower chamber pressures are reached. 
The first tangential (IT), second tangential (2T), and first radial (1R) modes are unstable over 
wide changes in operating conditions. An orifice diameter change from 0.09 to 0. 10 in. was 
evaluated to determine the effect of orifice diameter on stability. The increase in orifice diameter 
did not improve stability margin sufficiently to provide any significant advantages. Damping 
devices were predicted to be very effective in damping the unstable modes. Also, this design 
was very different from previously tested hardware in terms of number of elements and orifice 
size. 


Table 3 also provides a rating of the candidate injector designs. Using the criteria in Table 
3, the fine triplet injector was the best choice to validate the models contained within ROCCID. 
The fine triplet injector stability was sensitive to chamber pressure and mixture ratio variations 
that allowed for testing the sensitivity and predictive capabilities of the models. The fine triplet 
injector was sufficiently different from previously tested injectors to allow a true "a priori" 
prediction. As a result, this injector was selected for the validation test hardware. 

The test validation engine designed using the ROCCID methodology is shown installed on 
the Aerojet E-4 test stand in the photograph provided as Figure 4. A picture of the injector face 
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a. Predicted Stability Without an Acoustic Resonator 



Figure 3. Expected Combustion Stability For the ROCCID 
Validation Hardware 
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Figure 4. The ROCCID Test Validation Engine Assembly Installed on Test Stand E-4 


II. Executive Summary, (cont) 

is shown on Figure 5 and the OFO triplet element drawing is shown on Figure 6. The injector 
contained 105 OFO triplet elements with equal 0.090 in. diameter oxidizer and fuel orifices and a 
35 deg impingement half angle. 

Validation Test Results 

Validation tests were successfully conducted in the Spring of 1991 in the Aerojet E-4 test 
facility. A total of twenty-seven tests produced useable combustion stability and/or performance 
data. 


The actual test points are shown in Figure 7, identified by test number. The tests included 
eleven with a bituned acoustic cavity (Block 1 tests), eight with no cavity (Block 2 tests), and 
eight with a monotuned cavity (Block 3 tests). A list of individual test conditions, planned and 
actual, is provided on Table 4. As can be seen, the ranges of chamber pressure and mixture ratio 
tested correspond well to the originally intended ranges. The broad range of operating conditions 
provided an excellent test database for model assessment. 

The testing consisted of "3 blocks" of tests. Block 1 tests were with a bituned acoustic 
cavity and were structured to provide test data in the most combustion stable engine 
configuration. These tests established performance values and mapped the regions of stable and 
unstable combustion. They were predicted to be the safest tests from a hardware damage 
standpoint since high amplitude combustion instabilities were not likely as seen in Figure 3b. 

This provided the most benign environment to establish start and shutdown sequences and 
performance values. The Block 2 tests were without any acoustic cavity and provided baseline 
data on chamber sound speed and an direct assessment of the benefit of the cavity used in the 
Block 1 tests. The Block 3 group of tests provided test data with a monotuned acoustic cavity. 
The tests in these three blocks are shown on Table 4. 

An overview of the test logic for each of the Blocks of tests is shown on Figure 8. The 
path depicted by the shaded boxes was the logic path followed during the test program. The 
Block 1 tests with the bituned cavity were structured to meet the objectives of test series A (start, 
shutdown and flowrate balance verification), B (verification of stable operation at the nominal 
operating point, MR=2.8, Pc=1250 psia), C (verification of chug and 1L mode predictions), and 
D (verification of stable operation at nominal conditions). The Block 2 tests without acoustic 
cavities were structured to meet the objectives of test series E (verification of unstable operation 
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Figure 5. The Injector Core With Its Faceplate Installed for Brazing 
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Figure 6. The OFO Triplet Element Was Designed to ROCCID Specification 



LOX Flowrate, tonVsec 


Figure 7. ROCCID Validation Test Points 
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TABLE 4. PLANNED AND ACTUAL TEST CONDITIONS 



Planned Conditions 

Actual Conditions 



Chamber 


Chamber 


Test 


Pressure 

Mixture 

Pressure 

Mixture 

No. 

Block 

psia 

Ratio 

psia 

Ratio 

4 

I 

1250 

2.80 

1177 

2.59 

7 

I 

1500 

1.50 

1441 

1.45 

9 

I 

500 

2.80 

505 

2.94 

12 

I 

300 

2.80 

250 

1.01 

13 

I 

1250 

1.20 

1210 

1.13 

14 

I 

1500 

2.80 

1358 

2.67 

15 

I 

1000 

2.00 

969 

1.93 

17 

I 

800 

7.50 

792 

6.74 

18 

I 

1250 

2.80 

1165 

2.35 

19 

I 

1250 

2.20 

1220 

2.06 

20 

I 

1250 

1.50 

1193 

1.50 

21 

II 

1250 

2.80 

1222 

2.67 

22 

II 

1250 

2.80 

1092 

2.52 

23 

II 

1000 

1.20 

1004 

1.25 

24 

II 

800 

7.50 

788 

6.71 

25 

II 

1500 

2.80 

1410 

2.86 

26 

II 

1250 

1.50 

1130 

1.70 

27 

II 

1750 

2.80 

1706 

3.03 

28 

II 

1250 

5.00 

1260 

5.15 

29 

III 

1250 

2.80 

1208 

2.55 

30 

III 

800 

2.80 

791 

3.09 

31 

III 

1000 

2.00 

953 

2.09 

36 

III 

250 

2.80 

249 

3.31 

37 

III 

1750 

2.80 

1735 

3.05 

38 

III 

1500 

2.80 

1467 

3.07 

39 

m 

1250 

5.00 

1235 

5.60 

40 

hi 

1250 

2.20 

1174 

2.23 
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II. Executive Summary, (cont) 


without acoustic cavities at the nominal design point), F (verify stable operation without acoustic 
cavities at off nominal conditions, and G (verify unstable operation without acoustic cavities at 
off nominal conditions). The Block 3B tests with a reconfigured (monotuned) acoustic cavity 
were conducted since the original bitune cavity did not provide as large a region of stable oper- 
ation as originally desired.* This block of tests met objectives B (confirm nominal MR-Pc 
operation with a new cavity tune), D (confirm off nominal stability with a new cavity), and H 
(confirm off nominal unstable points with a new cavity). 

Acoustic Mode Stability 

The ROCCID methodology was used to recalculate the combustion stability for the tests 
where the acoustic cavity design was inadvertently changed (Block 1 and Block 3). Pretest 
predictions using ROCCID for the configuration without a cavity were still valid. The ROCCID 
calculations of combustion stability for each test condition are summarized in Table 5. Test 
results are also shown in Table 5 for comparison. 

The eleven Block 1 tests with a bituned acoustic cavity covered a range of chamber pres- 
sures from 250 to 1442 psia and mixture ratios from 0.54 to 6.74. Seven of these tests were 
stable and four had spontaneous first tangential mode instabilities that ranged from peak-to-peak 
amplitudes of 46 to 103% of the steady-state chamber pressure. Four of these stable tests had the 
combustion process disturbed by bombs to promote combustion instabilities. Bomb over 
pressures ranged from 5 to 67%. Two stable tests were low chamber pressure tests for chug 
stability evaluation and were not bombed and one test had a bomb intended but did not produce a 
measurable overpressure. The ROCCID methodology correctly predicted all of the four 
tangential mode unstable tests but only two of the seven stable tests. In some cases, ROCCID 
predicted instability in the 1L mode, but no instabilities in this mode were observed. 

The eight Block 2 tests without an acoustic cavity had chamber pressures from 788 to 
1706 psia and mixture ratios from 1.25 to 6.71. Only two of these tests were stable, the lowest 
and highest mixture ratio tests. Both these tests had bomb over pressures of 5 and 11%. Of the 
six unstable tests all were in the first tangential mode and all were correctly predicted to be 


* Subsequently, an error in the design of the acoustic cavity block was found that resulted in a significantly larger 
radial entrance to the cavity. The larger entrance resulted in an effectively shorter cavity and higher cavity gas 
temperature that considerably reduced damping near the first tangential frequency (See Figures 54, 55, and 56). 
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COMPARISON BETWEEN CALCULATED RESULTS AND TEST DATA 
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(1) Chamber pressure below which chug instability is predicted for given test mixture ratio. 

(2) AL is the growth coefficient calculated by ROCCID as illustrated In Figure 18. 

(3) Calculated results using ROCCID as described in Section IV H. 


II. Executive Summary, (cont) 


unstable by the ROCCID methodology. The methodology predicted one of the stable tests 
correctly (the low Pc chug test) but did not correctly predict the other. 

Block 3 testing consisted of 8 tests with a monotuned cavity that were conducted over a 
chamber pressure range from 249 to 1735 psia and a mixture ratio range from 2.09 to 5.60. Six 
tests were unstable as predicted by ROCCID and two tests were stable, only one of which was 
predicted stable by ROCCID. 

The combustion stability test results from all three test blocks are shown in Figure 9 as a 
function of measured chamber pressure and mixture ratio. An approximate correlation for the 
stable/unstable transition as a function of chamber pressure and mixture ratio is also provided 
and compared to the ROCCID prediction. From these results, no change in this transition 
boundary was discemable with the use of the improperly tuned acoustic cavities. Some 
additional damping was evident when these acoustic cavities were used, however, based on the 
reduction in the amplitude of the instability measured with acoustic cavities compared to that 
measured without cavities. Also, a reduced experimental growth coefficient was estimated from 
tests with the acoustic cavities compared to tests without cavities. Similar trends were predicted 
by ROCCID as noted in Table 5. 

A comparison of the experimentally observed stability transition boundary (from Figure 9) 
to that predicted by ROCCID (IT lower boundary from Figure 3a) is shown in Figure 10. From 
this comparison, the zone of unstable operation expected based on the ROCCID results was 
somewhat larger than the experimentally determined zone. Interestingly, however, a close corre- 
lation between the ROCCID prediction and the observed test results was achieved near the 
design mixture ratio of 2.8. This may indicate that the stability analysis models provide better 
predictions at mixture ratios near "nominal" values where a majority of the historical database 
was generated and upon which the models were developed and verified. 

Nonacoustic (Chug Stability) 

All tests were chug stable at the tested operating conditions. A few tests briefly 
encountered chug instabilities during the early part of the start-up transient where the chamber 
pressures were extremely low. These instabilities, however, disappeared quickly as the mean 
chamber pressure increased to its steady-state value. These tests results confirmed the chug 
predictions made using the ROCCID methodology. 
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open symbol : stable 

solid symbol: unstable 

number 1-20: test with bituned cavities 

number 21-28: test with blocked cavities 

number 29-40: test with a monotuned cavity 


Figure 9. ROCCID Validation Test Results Indicate Stable and Unstable 
Operating Regions 
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Figure 10. ROCCID Predicted a Somewhat Large Zone of Unstable 
Operation Than Was Observed Experimentally 


II. Executive Summary, (cont) 


Performance 

The test measured performance efficiency based on either characteristic velocity ( %C *) or 
specific impulse (energy release efficiency or ERE) was 5 to 7% lower than that predicted by the 
ROCCID analysis. A comparison of the predicted and measured performance in terms of the 
percent C* and the thrust based energy release efficiency (ERE) is shown on Table 6. The trends 
of both the measured and predicted performance with mixture ratio are shown on Figures 1 1 and 
12 for ERE and %C*, respectively. In addition to the measured data showing significantly lower 
performance than the predictions, the mixture ratio for minimum efficiency is somewhat 
different between the predicted and measured values. The predicted minimum efficiency is 
97.8% and occurs at a mixture ratio of 2.2, while the measured minimum (91%C*, 93% ERE) is 
at a mixture ratio of approximately 2.9. The ROCCID code predicts a vaporization efficiency of 
99.8% and a mixing efficiency of 98% at a mixture ratio of 2.2. 

The vaporization efficiency at near stoichiometric mixture ratios appears reasonable based 
on a comparison between the predicted and measured static pressure profiles in the combustion 
chamber. A typical pressure profile comparison at near nominal conditions is shown on Figure 
13. However, at very low or high mixture ratios the ROCCID code over predicts the vaporiza- 
tion rate, as can be seen from the comparison between the predicted and measured chamber static 
pressure profiles shown on Figure 14. This over prediction is probably a result of the generalized 
length correlation used for vaporization estimates that is based on data and calculations assuming 
typical stoichiometric combustion temperatures. The generalized length model in ROCCID may 
need to be modified to reflect the lower vaporization rates at off nominal conditions either by 
using a more mechanistic model, or by modifying the generalized length correlation by including 
a combustion temperature correction factor. 

Since the vaporization efficiency at near stoichiometric mixture ratio appears valid, the 
lower experimental efficiency must be the result of a lower mixing efficiency. The performance 
mixing efficiency was calculated by ROCCID from an input unielement E m of 0.77 obtained 
from previously conducted cold flow tests. A unielement to multielement correlation from 
DROPMIX was used by ROCCID to determine the multielement E m (.87) for the 105 element 
OFO injector. A better match of the prediction to the data is obtained, however, if a multi- 
element E m of .75 is used, which would be indicative of no interelement mixing improvement. 
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TABLE 6. MEASURED AND PREDICTED VALIDATION ENGINE PERFORMANCE 



Calculated 

Measured 

Test 

Pc 

Mr 

Eta-C* 

ERE 

Eta-C* 

ERE 

4 

1178 

2.59 

0.983 

0.983 

0.801 

0.842 

7 

1441 

1.45 

0.990 

0.989 

0.975 

0.972 

9 

505 

2.94 

0.988 

0.988 

0.994 

0.987 

12 

250 

1.01 

0.981 

0.980 

0.863 

0.792 

13 

1210 

1.13 

1.003 

1.004 

1.049 

1.010 

14 

1460 

2.67 

0.985 

0.985 

0.932 

0.948 

15 

969 

1.93 

0.980 

0.981 

0.968 

0.951 

17 

795 

6.74 

0.991 

0.991 

1.025 

1.015 

18 

1165 

2.35 

0.979 

0.979 

0.978 

0.997 

19 

1220 

2.06 

0.979 

0.980 

0.994 

1.000 

20 

1193 

1.50 

0.983 

0.982 

0.955 

0.948 

21 

1222 

2.67 

0.984 

0.984 

0.934 

0.928 

22 

1092 

2.75 

0.980 

0.980 

0.836 

0.924 

23 

1004 

1.25 

1.014 

1.016 

0.961 

0.937 

24 

788 

6.71 

0.990 

0.990 

0.956 

0.947 

25 

1410 

2.86 

0.987 

0.987 

0.908 

0.926 

26 

1130 

1.70 

0.983 

0.983 

0.736 

0.753 

27 

1755 

3.03 

0.990 

0.990 

0.890 

0.938 

28 

1260 

5.15 

0.994 

0.994 

0.904 

0.946 

29 

1208 

2.55 

0.982 

0.982 

0.927 

0.915 

30 

791 

3.09 

0.991 

0.990 

0.919 

0.921 

31 

953 

2.09 

0.978 

0.979 

0.843 

0.663 

36 

249 

3.31 

0.990 

0.990 

0.931 

0.899 

37 

1738 

3.05 

0.990 

0.990 

0.918 

0.927 

38 

1467 

3.07 

0.991 

0.990 

0.920 

0.963 

39 

1214 

5.60 

0.994 

0.994 

0.945 

0.916 

40 

1180 

2.23 

0.978 

0.978 

0.874 

0.886 


Note: Efficiencies greater than 1 .00 are suspect and are probably the result of interpolation 

accuracies for calculated values or test measurement accuracies for measured values. 
See Section IV H. 
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Figure 13. Good Agreement Observed Between Predicted and Measured Chamber 
Static Pressure Profiles at Near Stoichiometric Mixture Ratio 
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Figure 14. Poor Agreement Observed Between Predicted and Measured Chamber 
Static Pressure Profiles at Very Low and High Mixture Ratios 


II. Executive Summary, (cont) 


Cold flow testing, shown in the photograph of Figure 15, indicated that the propellant 
spray fan formation from individual elements is not complete before adjacent elements begin to 
restrict the fans* . This may cause a multielement mixing efficiency to be actually lower than the 
unielement mixing efficiency. To confirm this, it is recommended that both unielement and 
multielement cold flow tests be conducted to measure the cold flow mixing efficiencies. The 
correlation currently in the DROPMIX code indicates an increase in multielement Em for a given 
unielement Em as the pattern fineness increases. This correlation should be used with caution for 
very fine patterns and an expanded element mixing database covering a wide range of element 
pattern densities and element types using both cold flow and hot fire data should be generated. 

At a minimum, both unielement and multielement Em's should be determined from cold flow 
testing for designs that are outside of the DROPMIX database. 

The validation test program has shown that ROCCID is a good tool that can be effectively 
used to perform first cut combustor design and parametric analyses. The testing showed some 
limitations within the methodology such as uncertain combustion response calculations, a 
constrained interelement mixing correlation and an apparent over prediction of the propellant 
vaporization rate at very off-nominal mixture ratio. However, the ROCCID stability predictions 
were conservative in that ROCCID tended to predict less stable operation than what was obtained 
and the overprediction in performance efficiency probably could have been found through 
relatively simple multielement cold flow mixing tests. An overall assessment of the ROCCID 
methodology and the results achieved during the validation test program are provided in the 
following section, III Conclusions and Recommendations. 


* Cold flow visualization tests of OFO triplet elements as a function of element spacing have been performed by the 
NASA Lewis Research Center and results supporting this observation are presented in NASA Technical 
Memorandum 105750 (AIAA-92-3226). 
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Figure 15. ROCCID Injector Cold Flow Propellant Spray Fan Formation 
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in. CONCLUSIONS AND RECOMMENDATIONS 


A computerized design methodology has been developed for use in designing efficient and 
stable combustors for LOX/Hydrocarbon and LOX/Hydrogen rocket engines. This Combustor 
Analytical Interactive Design Methodology (ROCCID) is unique in several respects. First, it 
contains an ordered assemblage of industry-wide analysis models that are used to define the 
effects of design and operating parameters on combustor performance and stability. Second, all 
practical and available analysis models and correlations are included even when redundant capa- 
bilities exist. This allows the user to directly compare the results of available techniques and 
select the method that appears to be the most applicable for a specific application. It also pro- 
vides a means to identify design and operational uncertainties based on the differences predicted 
by available techniques. Third, these analysis models and correlations are contained in a mod- 
ular format so that model updates and new models can be incorporated into the methodology as 
they are developed. Last, ROCCID is designed to be user friendly and contains an interactive 
front end for input formulation and checking, file generation, and output display. 

Validation tests using hardware with critical design features determined from the ROCCID 
methodology has demonstrated both the utility and capability of the methodology and identified 
some of its shortcomings. ROCCID high frequency (acoustic coupled) instability predictions 
proved to be conservative in that the experimental zone of instability was somewhat smaller, in 
terms of operating chamber pressure and mixture ratio, than that predicted from the ROCCID 
methodology. Nonetheless, the general unstable character of the combustion chamber without 
acoustic damping was adequately predicted from a design selection standpoint. In this case, use 
of the ROCCID methodology would have led to improved stability margin through the incor- 
poration of acoustic damping devices to increase chamber damping and/or a change in the 
injector design to reduce the predicted combustion response. 

In spite of the relatively good correlation between the ROCCID predictions and the valida- 
tion test results, determination of the combustion response remains a highly empirical and 
analyst dependent process. The adequacy of a-priori specification of the combustion response 
can be a function of the element type and design, the propellants and the available experimental 
database. While some analytical approaches, such as CRP and the NASA Lewis Injection 
Coupling Model have been developed, their application has been relatively limited and generally 
involved in correlation of tested designs. Clearly, further advancement in the modeling and 
characterization of the combustion response is essential for improvement in the overall predictive 
capability of liquid propellant combustion stability. 
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IE. Conclusions and Recommendations, (cont) 


However, since the whole sensitive time lag methodology, upon which ROCCID is based, 
is built upon empirical correlations and analytical approaches that have been verified, refined 
and/or anchored with empirical data generated over the past thirty to thirty-five years, it is not 
obvious that additional improvements in any specific modeling area would result in an overall 
improvement in the combustion stability characterization. For example, the chamber response 
modeling does not account for the acoustic energy damping provided by the chamber wall and 
the propellant spray. Improvements in this area may be able to be incorporated into the chamber 
response model to make it more physically representative. Yet since, the combustion response, 
n, has been "backed out" from experimental results using chamber response models without the 
wall and spray damping considerations, its use with the improved chamber response models may 
actually result in a degradation in the predictative capability of the methodology. 

Chug stability limits appeared to be adequately defined. The predicted longitudinal insta- 
bility at lower chamber pressure was not observed at steady-state conditions but was observed 
during the start transient. Thus injection pressure drop requirements can be adequately defined 
using the ROCCID methodology. 

Performance efficiency was overpredicted by ROCCID. In this case, it appears that the 
interelement mixing efficiency was overpredicted by a basically empirical correlation. The 
propellant vaporization appeared to be adequately predicted at near nominal mixture ratios but 
was overpredicted at very off-nominal mixture ratios where the driving temperature (combustion 
temperature) is significantly reduced. In this sense, the validation test result can be used to 
extend the empirical limits of the ROCCID models for very off-nominal mixture ratio operation 
with LOX/RP- 1 propellants. Additionally, future improvements may result from emerging 
computational fluid dynamic (CFD) modelling that can, as a minimum, be used to validate or 
improve existing and primarily empirical based techniques while using ROCCID as a parametric 
design tool to home in on specific areas to study with CFD. 

ROCCID will provide significant cost savings for the design and development of rocket 
engine combustion chambers. Standardization of the design and analysis approach and the avail- 
ability and use of all the best combustion analysis models and correlations ensure cost efficiency 
improvements. ROCCID is a low cost design tool in that it results in relatively low CPU usage, 
doesn’t require high performance computers such as a CRAY, is portable to different computer 
platforms and provides, through its built-in model connectivity and interactive front-end, 
increased parametric design and analysis capability with decreased engineering hours. Savings 
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III. Conclusions and Recommendations, (cont) 


of 50% or greater in manhours are estimated when using ROCCID once the analyst is familiar 
with its operation and capabilities. 

Further validation and continued wide application experience will provide additional 
means to identify strengths and limitations in the methodology. Isolation of weaknesses will 
provide for focused improvements in the methodology and an evolution to a proven and robust 
combustor analytical design technology that can serve the industry while sophisticated and 
mechanistically based computational methods are developed. 
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IV. TECHNICAL NARRATIVE 


The LOX/Hydrocarbon Rocket Engine Analytical Design Methodology Development and 
Validation was a three year program conducted by the Aerojet Propulsion Division under con- 
tract NAS 3-25556 from the NASA Lewis Research Center. This program consisted of two 
major efforts: first, the development of an analytical design methodology for liquid propellant 
rocket engines that uses existing industry analysis computer programs and empirical correlations, 
and second, the validation and critical evaluation of this methodology through its application for 
the design of a test demonstration thrust chamber and prediction of its combustion stability and 
performance characteristics. The capability of this methodology was demonstrated through fab- 
rication and testing of this thrust chamber over a wide range of mixture ratio and chamber pres- 
sure while obtaining detailed measurements of its operation for comparison to the performance 
and combustion stability forecasts developed using the methodology. 

These activities were completed during eight technical tasks, the results of which are sum- 
marized in this final report. Task 1 .0 consisted of a survey of existing industry analysis models 
for use in the methodology and existing stability and performance databases for evaluation of 
these models. Development of the methodology, which was named the ROC ket Combustor 
Interactive Design (ROCCID) methodology, was completed during Task 2.0 and a limited 
release of a test (beta) version of the computer code was made. The newly developed ROCCID 
methodology was applied to illustrate the effects of critical design parameters and operating 
conditions on combustion chamber performance and combustion stability during Task 3.0 to 
provide a basis for the selection of these parameters for hardware design to validate the 
methodology. Task 4.0 consisted of the development of a detailed test plan and thrust chamber 
performance and combustion stability predictions for validation testing of the newly developed 
ROCCID methodology. Validation hardware design and hardware fabrication were completed 
during Tasks 5.0 and 6.0, respectively. Finally, hot fire testing of the validation hardware was 
completed during Task 7.0, and data analysis and comparison of test results with the ROCCID 
predictions was provided in Task 8.0. A summary of these task results are contained in the 
following subsections. 

A. DATA/STABILITY MODEL SURVEY 

A survey of existing analysis models and available databases containing combustion 
stability and performance measurements was conducted to identify and acquire analysis models 
that perform critical analytical characterization or empirical correlation of important combustion 
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IV.A, Data/Stability Model Survey, (cont) 

processes that influence the combustion stability and performance of a rocket engine combustor. 
Models reviewed covered three general areas of combustion analysis; (1) steady-state combus- 
tion, which models the propellant injection, atomization, vaporization and mixing processes, (2) 
combustion stability, which models the non-steady injection, combustion and chamber response 
of a combustor design and (3) performance, which defines the combustion efficiency of the 
energy release processes. 

Evaluation criteria were developed to provide a means for rating and selecting 
analysis models for further review and comparison with the selected combustion databases. The 
evaluation criteria used to rate various analysis models are provided in Table 7. The first three 
criteria were considered most important since a failure to meet any of these three requirements 
would mean that a particular model was unacceptable for use in the analytical design methodol- 
ogy. For this reason, the relative value of these criteria range from 0 to 1.0 and these are used as 
multipliers in the overall rating value equation also shown in Table 7. A zero value for the A, B, 
and C evaluation criteria will result in a zero value for the overall rating value. In this case it is 
clear that a model is unacceptable for use in the advanced methodology. Conversely, a value of 
1.0 for all three of these criteria means that selection of the model for further evaluation will be 
made on the basis of the remaining three evaluation criteria. 

The remaining three evaluation criteria, D, E, and F, have relative values between 1 
and 10. The higher the relative value for a particular model, the better that model meets the 
evaluation criteria. Note that the average relative value of the D, E, and F evaluation criteria are 
used in defining the overall rating. 

A list of 18 models that were evaluated for further study is contained in Table 8. 
These models are grouped in three categories: (1) Steady State Combustion, (2) Combustion 
Stability, and (3) Engine Performance. The table provides our subjective assessment of the rela- 
tive value for each of the six evaluation criteria for each model. An overall assessment is 
included using the overall rating value equation from Table 7. The final column in Table 8 indi- 
cates if a model is recommended for further evaluation by correlation with existing industry 
databases. A favorable recommendation is made if the overall rating value is 5.0 or greater. 

The set of recommended models included in Table 8 is sufficient to build a complete 
analytical design methodology with only one exception. The ability to model combustion insta- 
bility coupling with nonlinear feed system acoustics is not contained in the set of recommended 
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TABLE 7. ANALYSIS MODEL EVALUATION CRITERIA HAVE BEEN DEVELOPED 
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Overall Rating Value = (A)(B)(C) 


TABLE 8. COMPUTER CODE EVALUATION SUMMARY 
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IV.A, Data/Stability Model Survey, (cont) 


models. This deficiency was eliminated through the development of non-linear injection model 
(Leinj) similar to the approach described in Reference 2. 

Existing combustion performance and stability databases listed in Table 9 were also 
reviewed as a source for data and design information for further model evaluation. Four of the 
eleven data bases reviewed were selected for this evaluation.* These sets cover a range of hydro- 
carbon fuels (and hydrogen) injector elements (impinging and concentric tube) acoustic damping 
devices and operating conditions (Pc, O/F, Dchb, Tfuel, etc.). Most importantly, they contain 
stable and unstable operating characteristics and performance data that can be compared to 
model calculations. 

Unfortunately, not all of the 1 1 recommended models were evaluated and compared 
to the existing data bases. Model availability, computational flaws, inefficient run times, ill- 
defined input and extensive and previously documented model evaluations were reasons for 
exclusion of some models. Nonetheless, a complete set of models were available to calculate 
combustor performance and combustion stability. In any case, several critical model correlations 
such as spray dropsize or burning response were available and were compared to each other and 
to data base measurements. 

Performance model evaluations were made by comparing calculated combustor effi- 
ciency (energy release or C*) to test based published values. Calculated variations in efficiency 
with key combustor design parameters such as injector element configuration and chamber 
length and operating parameters such as fuel type and mixture ratio were compared with experi- 
mentally determined trends, combustion stability models were evaluated by comparing predicted 
stable or unstable operation with documented results from the selected data bases. Chug, 
burning coupled and injection coupled instabilities were modeled and acoustic mode and 
frequency predicted. Injection element type, chamber size, acoustic damping devices (resonators 
and baffles) and fuel type are parameters that were included in this evaluation. 

Results from the performance model evaluation have shown that combustor effi- 
ciency can be predicted to within 1 to 2% over a wide range of design and operating conditions. 
Injector designs producing high efficiency, such as the impinging OFO triplet element, were well 


* Subsequently, data anomalies have been discovered in the LeRC in-house data base (NASA TM 79319) and 
further use of this data base is not recommended pending its correction. 
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TABLE 9. DATA BASE REVIEW AND SELECTION 
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IV.A, Data/Stability Model Survey, (corn) 


characterized. Moderate to low performance designs, such as the like doublet with little or no 
intentional unlike propellant fan impingement, were not as well characterized because of defi- 
ciencies in modeling their propellant mixing inefficiencies. The injection process of gas/liquid 
concentric tube element injectors was found to be not well characterized. However, because of 
the volatility of liquid oxygen (the fuel is injected at a temperature above its critical value with 
this injection concept) the vaporization lag produced with concentric tube injector is relatively 
small and therefore less sensitive to LOX atomization model assumptions. Similarly, the intra- 
element mixing process of this element is not well characterized but mixing inefficiencies have 
been overcome by maintaining small mixing distances through the use of relatively small ele- 
ments. Performance model limitations could become more significant for this element type with 
the use of higher molecular weight fuels (e.g., methane vs hydrogen) and larger (fewer and lower 
cost) injection elements. 

The propellant vaporization rate or the droplet heating rate (for liquids above their 
critical pressure) appears to be adequately characterized using the simplified "generalized length" 
correlation developed by Priem in NASA TR-67. This is particularly true in light of the uncer- 
tainties in the propellant atomization models. Use of generalized length correlation will greatly 
streamline the propellant vaporization calculations and improve the efficiency of the 
methodology. 

Stability model evaluations compare surprisingly well with the selected data base 
results. This is particularly true for the liquid/liquid impinging element injector designs with and 
without acoustic damping devices. The stabilizing effects of both baffles and acoustic cavities 
were well characterized when compared with test results. Burning rate coupled instabilities were 
predicted with likely acoustic modes and frequencies corresponding closely with experimental 
results for different chamber sizes, element designs and orifice sizes. The burning response was 
found to contribute the largest uncertainty in the model predictions of this mode instability. 
Burning response, based on the empirical pressure index (n) and sensitive time lag (t) 
correlations, has been previously characterized as accurate to within only ±25%. Nevertheless, 
when used with stability model chamber admittance predictions, good insight into combustion 
stability of individual designs was achieved. An analytically based combustion response model 
(CRT) also appears promising although its application in this evaluation was somewhat limited. 
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IV.A, Data/Stability Model Survey, (corn) 


The combustion stability of the LOX/CH 4 shear coaxial element combustor was not 
well characterized using available models. Model characterization of the combustion stability 
was improved by using the experimental results to infer or anchor the model(s) of the element 
injection process. In this case, the combustion stability could be shown to be injection rate 
coupled rather than burning rate coupled and therefore a strong function of the fuel injection 
velocity (mixture ratio and methane density). Improvement in the concentric tube injection 
process modeling is essential before good a priori predictions of its combustion stability can be 
made. Details of the performance and stability model correlations are provided in Reference 3. 

Finally, the models and correlations evaluated were reassessed on the basis of the 
model evaluation results and a complete set of models was recommended for incorporation into 
the rocket combustor analytical design methodology. In many cases, more than one model or 
correlation was recommended to provide the same capability or function so the user has options 
for model selection and may produce comparative solutions (predictions) upon which to base his 
design decisions and assessment of critical design and operating parameters. A listing of the 
models included in ROCCID is provided in the following subsection. 

B . THE ROCKET COMBUSTOR INTERACTIVE DESIGN METHODOLOGY 

The design and/or evaluation of a liquid propellant rocket engine combustion 
chamber is normally achieved through an interactive process consisting of inter-dependent 
design analyses. These design analyses cover a multitude of disciplines including aerodynamics, 
thermodynamics, combustion, hydraulics, heat transfer and structural analysis. Many 
subcomponents are involved in creating a combustion chamber design including an injector and 
associated propellant manifolds, a combustion chamber and its associated cooling schemes, a 
supersonic expansion nozzle, an ignition system for non-hypergolic propellants and possibly 
acoustic damping devices such as resonators and/or baffles for combustion stability control. 

A flow chart containing a description of the first order analyses and design iteration 
for the combustion design process is shown in Figure 16. A portion of this methodology, noted 
with a dashed line boundary in Figure 1 6, has been developed into a user friendly interactive 
computerized methodology. This methodology, entitled ROCket Combustor Interactive Design 
or ROCCID, contains the design and analysis logic and computer analysis modules to define the 
combustion stability and steady state performance of a combustor. 
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Figure 16. ROCCID Performs a Majority of the Combustor Design Iterations 
to Meet Combustion Stability and Performance Requirements 
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IV.B, The Rocket Combustor Interactive Design Methodology, (cont) 

The structure of ROCCID is illustrated in Figure 17. ROCCID contains three main 
components including: 

(1) An interactive front end (IFE) that assists the user with input setup, generation 
of plots of model output data, maintenance of combustion gas properties and 
logic for controlling file structure for the reuse and restart of computer runs. 

(2) A point analysis option that provides performance and combustion stability 
analysis of existing combustor designs. 

(3) A point design option that creates essential design features of a combustor that 
satisfies specified performance and combustion stability requirements. 

Both the point analysis and point design options use a library of performance and 
combustion analysis models that were selected from existing appropriate codes. These analysis 
models are contained within ROCCID in a modular format. This permits the user to access 
specific models for a specialized subanalysis or to use two or more models that perform similar 
functions to define and resolve uncertainties in that particular area of the analysis. Modular con- 
struction also permits easier and less costly upgrades to the methodology as new analysis models 
are developed. 

The IFE is a menu driven preprocessor constructed using an extensive library of 
interactive subroutines. Each input character is checked in the IFE for validity. Warnings are 
displayed when input errors are sensed. Replay files that contain all case input are created and 
maintained. These files can be edited and used as input for a subsequent session. The user may 
repeatedly alter the input until the desired result is achieved. Required propellant properties for 
both the injected fuel and oxidizer and their resulting combustion gas are internally generated 
and maintained in files. The "MIPROPS” data base, generated by the National Bureau of 
Standards (NBS), is used to obtain properties for LOX, H2, CH4 and C3H8. RP-1, N2O4, and 
MMH properties have also been included using special curve fit correlations. Finally, 
combustion gas properties are obtained using the One Dimensional Equilibrium (ODE) code 
developed by the NASA Lewis Research Center. 

Upon completion of the input deck, the user may execute either the point analysis or 
point design options. These options contain many interactive decision points for the user. Upon 
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Figure 17 . Major Components of ROCCID include the Interactive Front End, 
Point Analysis, and Point Design Evaluation 
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IV.B, The Rocket Combustor Interactive Design Methodology, (cont) 

completion of the analysis or user termination, the run stream returns to the IFE to monitor the 
results, record, print pertinent information, and prepare input for the next run. 

The performance of the combustor is defined by the energy release efficiency. This 
accounts for combustion efficiency limitations resulting from incomplete propellant vaporization 
and/or mixing. The energy release efficiency is calculated using the JANNAF simplified per- 
formance procedure (Ref 4) and the propellant vaporization and mixing limitations from the 
steady-state combustion analysis. An input file for use in the TDK/BLM performance analysis 
of the complete engine is also generated. 

A steady-state combustion analysis including propellant atomization, vaporization, 
and mixing supplies key input to the performance and stability analyses. Four 
models/correlations for propellant dropsize are included for standard injector elements including 
showerhead, doublet and triplet impinging elements, and shear and swirl coaxial elements. Drop 
sizes from all applicable correlations are calculated and displayed for comparison. A user may 
select any of the calculated values for performance and stability analyses or provide another 
estimated value. Propellant vaporization is calculated using the generalized length correlation 
(Ref 5). Propellant mixing is based on the use of a unielement mixing efficiency value deter- 
mined previously from cold flow measurements and adjusted for interelement mixing effects. 

The combustion stability analysis can be made with a large array of models used to 
calculate the chamber admittance and the burning and injector responses magnitudes. These 
models provide the capability to estimate combustion stability margin for all common types of 
combustor instabilities including chugging and chamber acoustic coupled (high frequency) 
modes. The effects of damping devices such as acoustic cavities and baffles are also considered 
by these models. A summary of the combustion stability models available within the ROCCID 
is provided in Table 10. 

The chug stability characteristics are determined by calculating an acoustic chamber 
response without damping devices. This response is compared to a lumped parameter or acoustic 
injection response (user option) to determine the marginal chamber pressure for the on-set of 
chug stability. The injection response is determined as a function of operating pressure by 
evaluating combustor operation through flowrate reduction at a constant mixture ratio. 
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IV.B, The Rocket Combustor Interactive Design Methodology, (cont) 

High frequency stability is calculated using the sensitive time lag approach. 
Analysis models are available to calculate the chamber response with or without acoustic 
damping devices and both burning and injection responses. Table 10 identifies the three models 
currently available to estimate the chamber response. The HIFI model (Ref. 6) developed by 
Aerojet calculates the chamber response for chambers without acoustic damping devices or with 
monotuned or bituned 1/4 wave tube or helmholtz acoustic resonators. A second model, 

DIST3D (Ref. 7), developed by Colorado State University is also available. This model 
calculates the chamber response including damping generated from radially oriented baffles. 
Acoustic resonators are also modeled as secondary damping devices. A third model, FDORC 
(Ref. 8), is also available. This model calculates chamber response for piecewise distributed 
combustion with arbitrarily located 1/4 wave tube or Helmholtz resonators. 

Multiple approaches are also available to calculate both the burning response 
and the injection response as noted in Table 10. Pressure interaction index sensitive time lag 
(N/t) correlations developed by Reardon are included for the impinging doublet and triplet ele- 
ments and the concentric tube (coax) element. In addition, a combustion response model (Ref. 9) 
recently developed by Aerojet is included. This model is based on the Agosta-Hammer non-lin- 
ear vaporization response model. Two models are also available to calculate the injection 
response of the system. The first model (INJ) bases the response on element injection timelags 
including atomization and partial vaporization. The second model (developed by Priem and 
Breisacher) is appropriate for shear concentric tube elements and includes the effects of acoustic 
waves in the fuel and LOX injection manifold (or L OX tube) in the injection response. 

ROCCID displays the stability analysis results in terms of the calculated growth 
coefficient (X) for the particular acoustic mode of concern. This growth coefficient represents 
the amount of amplification required by the chamber to achieve the condition where the driving 
required to support the waveform exactly equals the driving response present within the system. 
Thus a value of X = 0 represents a neutral stability condition, X > 0 represents an unstable condi- 
tion and X < 0 represents a stable operating condition. 

The growth coefficient, in terms of the traditional combustion instability transfer 
function format, is illustrated in Figure 18. Three typical situations are illustrated in this 
example. Figure 18 (a) illustrates results of a traditional stability analysis where the region 
formed by overlapped chamber response and combustion response curves indicate the potential 
for an instability. However, this traditional analysis does not include a survey of chamber 


RPT/F0213.105-IV 


51 


9/30/92 


Magnitude 


Chamber 

response, 

X 


0 

500 

600 

-600 

Combustion response 



/ 

Region of potential 
unstable operation 


(a) Traditional neutral 
stability analysis, not 
quantitative. 




(c) Stable (X < 0). 


* 


# 


Figure 18. Illustration of the Growth Coefficient, X, in Terms of Traditional 
Combustion Stability Transfer Function Format 
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IV. B, The Rocket Combustor Interactive Design Methodology, (cont) 

admittances that would lead to a growth coefficient value that would provide zero stability 
margin, i.e., chamber response and combustion response curves being tangent at some point. 
The growth coefficient at this zero margin point provides a quantitative assessement stability. 
Figures 18 (b) and (c) illustrate unstable and stable quantitative predictions, respectively, where 
the growth coefficient is used to define the relative stability based on the rate of pressure 
amplitude growth (X. > 0) or decay (A,<0). The following is the exponential relationship used to 
define X: 


A/A 0 = e^ 1 

where A/Ao is the pressure amplitude ratio chamber over the period. At (sec). 

The resultant growth coefficient can be related to the damp time required for 
dynamic stability as defined in CPIA 247, by defining a growth coefficient (Acpia 247) required 
to achieve a 10:1 overpressure damping within the damp time specified from CPIA 247. That is: 

^CPIA 247 = In (0.1)/tdamp> CPIA 247 


where 

tdamp, CPIA 247 = 1.250/Vf (1/sec), 


and 


f = frequency of the unstable oscillation 
Reducing results in the following relationship: 

^CPIA 247 = -1.842 Vf (1/sec) 

Note that Acpia 247 will always be negative, indicating a damped system, and will be 
a function of the frequency (mode) of the damped oscillation. Therefore, if the growth 
coefficient is equal to or less than Acpia 247 the system is dynamically stable. Conversely, if the 
growth coefficient is greater than XcPIA 247 the system is not dynamically stable. 

The design requirements for combustor cooling are established using techniques 
outside of ROCCID. These requirements may include estimates of fuel film cooling required for 
chamber and baffle walls, dump cooling of baffle tips, and bulk temperature increases resulting 
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IV.B, The Rocket Combustor Interactive Design Methodology, (cont) 

from regenerative cooling of the nozzle, chamber, and resonator/baffle components. This infor- 
mation is used to calculate the propellant injection temperatures, injection orifice requirements, 
and the local flow injection mixture ratios. 

Capabilities and Limitations 

ROCCID has been specifically formulated to be applicable to combustor designs for 
LOX/HC and LOX/H2 propellants. Propellant and combustion gas properties for LOX-RP-1, 
LOX/CH4, LOX/C3H8, LOX/H2, and N2O4/MMH are internally calculated. Both gas/liquid and 
liquid/liquid propellant injection are considered. Conventional impinging like doublet and triplet 
(OFO and FOF) elements, unlike doublet elements, non-impinging showerhead and shear and 
swirl coaxial elements are modelled. The injector can consist of a mixed element pattern, 
including core, baffle, barrier and fuel film/cavity cooling elements. Different element types can 
exist in different zones (i.e., baffle core, barrier). However, in any one zone only one element 
type may exist and doublet elements must be in matched pairs of oxidizer and fuel. 

The point design option permits the user to constrain some design parameters, such 
as contraction ratio (CR) and chamber length (L'), thereby focusing on the best trade-off between 
performance and combustion stability. A simple trade-off between nozzle length and chamber 
length is also included to optimize engine delivered specific impulse for an envelope limited 
system. Acoustic damping devices are also recommended and their design features specified to 
provide the required combustion stability margin. Design trade-offs for a throttling engine are 
also performed. 

The point design option features an optimization of the injector element design. The 
guidelines in Table 1 1 are provided to aid the user in injector element type selection for a par- 
ticular application. The quantity of elements and the injector orifice size are calculated through a 
series of trade studies to satisfy performance and high frequency combustion stability require- 
ments, chug stability and other user-specified constraints. 

ROCCID has been prepared with certain limitations in order to simplify its construc- 
tion and guarantee is timely and affordable initial development. No supersonic nozzle effects are 
currently included. The nozzle design and engine specific impulse must be determined outside 
of ROCCID. Precombusted fuel (staged combustion cycle) is not presently considered. Mass 
addition from ablation, igniters, gas generator dump or transpiration cooling is not modelled. As 
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TABLE 11. Element Characteristics for Several Conventional Liquid Rocket Injector Elements Modeled Within ROCCID 
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IV.B, The Rocket Combustor Interactive Design Methodology, (cont) 


previously noted, cooling requirements are defined outside ROCCID, but their effects on per- 
formance and stability are considered. 

ROCCID has been developed and is operational on VAX 8600 series computers at 
both Aerojet and NASA Lewis Research Center. It has also been operated on a SUN4/SPARC 
architecture computer at SEA, Inc. Interactive graphics for Tektronix 40XX and 41XX terminals 
are also provided. The code has been constructed without machine dependent instructions, but 
operation on other computer systems remains to be verified. A complete description of the 
ROCCID methodology and a user's manual for its operation is provided in References 10 and 11. 

C. APPLICATION OF THE ROCCID METHODOLOGY 

The creation of a combustor design consists of an evolutionary process where design 
requirements and operational goals are used to establish design parameters for the injector and 
combustion chamber. Iterations to the set of critical design parameters are performed to meet 
performance goals while providing for stable combustion. ROCCID has been used to examine 
the sensitivity and influence of several injector and combustor design parameters on the 
predicted combustor performance and combustion. 

The propellant combination of LOX/RP-1 was selected for study using ROCCID and 
for eventual hot fire testing to validate ROCCID capabilities. This propellant was selected 
because historically, it has proven to be a difficult propellant to provide both stable combustion 
and high performance. 

Nominal design propellant injection temperatures of 174.4°R for liquid oxygen and 
ambient conditions (515°R) for RP-1 were selected. These conditions represent normal 
propellant storage temperatures within the Aerojet test facility and therefore were the most cost 
effective test conditions. 

A nominal design mixture ratio of 2.8 was selected since this approximates an 
optimum performance design from a thrust chamber performance viewpoint. Note that this 
mixture ratio value was selected for design purposes only. The effects of mixture ratio variation 
on a fixed design were evaluated using ROCCID and mixture ratio was a key operating 
parameter for validation testing. 
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IV.C, Application of the ROCCID Methodology, (cont) 

Thrust chamber size was also a key design parameter, particularly from a combustion 
stability standpoint. The chamber diameter has a direct effect on the resonant frequency of 
acoustic modes within the chamber, large chambers appropriate for booster engine 
application(Dc = 17.5 to 44 in.), have lower resonant frequencies that are more likely to couple 
with typical combustion responses. On the other hand, subscale chamber diameters of 5.5 and 
7.68 in. have been successfully used in the past for combustion stability investigations (Refs. 12 
and 13). For this study, a nominal baseline chamber diameter of 7.68 in. was selected since it is a 
proven "subscale" design and residual hardware was available from the test program in Ref. 13. 

A contraction ratio of 2:1 that corresponds to a throat diameter of 5.43 in. and a 
nominal design chamber pressure of 1250 psi were selected to provide a large operating chamber 
pressure variation within the existing hardware and test facility capabilities. The operational 
envelope of the validation hardware for the Aerojet E-4 test facility is shown in Figure 19. A 
variation in mixture ratio from approximately 1 to 10 was possible at the nominal chamber 
pressure of 1250 psi. A maximum chamber pressure of approximately 1800 psi could be 
achieved at the nominal mixture ratio. The mixture ratio excursion range diminishes as chamber 
pressure is increased due to the facility run tank pressure limits. 

Impinging injector designs were selected for this evaluation since they are most 
appropriate for liquid/liquid injection. Specifically, the like doublet and OFO triplet elements, 
both of which are included in the ROCCID methodology were evaluated. The OFO triplet 
element was selected over the FOF since the propellant density and mixture ratio resulted in 
equal orifice diameters with an OFO element but not with the FOF triplet. From experience, 
equal orifice diameters are desirable with a triplet element to provide optimum mixing and 
atomization. 

A summary of the selected design parameters for this study is presented in Table 12. 
Design parameters from two recent combustion stability investigations in Refs. 12 and 13, are 
also provided for comparison purposes. 

Downselect Criteria 

Criteria were established to select the injector for use in the validation of ROCCID. 
Because of limited resources, one and at best two different thrust chamber designs could be 
selected for design and testing. The following criteria were established on the basis of providing 
the best test of the capabilities of ROCCID: 
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Figure 19. Facility Operating Limitations 
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TABLE 12. COMPARISON OF SELECTED NOMINAL DESIGN PARAMETERS 



Combustion Stability Investigation 


NAS 3-255560) 

F0461 1-85-C-OlOOl 1 2 3 ) 

NAS 3-24612(3) 

Propellant Combination 

LOX/RP-1 

LOX/RP-1 

LOX/CH4 

Mixture Ratio 

2.8 ± TBD 

2.9 ± 0.3 

-1.6 

3 - 5 +0.2 

Injection Temperature 
LOX °R 
Fuel °R 

174 

515 

174 

530 ±20 

174 

530 to 438 

Chamber Pressure (psia) 

1250 ± TBD 

+100 
2000 _ 130 o 

+100 
2000 -500 

Chamber Diameter (in.) 

7.68 

7.68 

5.66 

Throat Diameter (in.) 

5.43 

4.43 

3.31 

Contraction Ratio 

2:1 

3:1 

2.92:1 

Chamber Length (in.) 

TBD 

8.5, 13,20 

14 

Nominal Thrust Level 
(Ibf) 

52,000 

54,000 

30,000 

Total Flowrate (lbm/sec) 

160 ± TBD 

+10 
180 -90 

+7 
90 -15 

Injection Element Type 

OFO Triplet 
Like Doublet 

OFO Triplet 
Like Doublet 

Shear Coaxial 


(1) LOX/Hydrocarbon Rocket Engine Analytical Design Methodology Development and 
Validation 

(2) Oxygen/Hydrocarbon Injector Characterization (Ref 13) 

(3) LOX/Hydrocarbon Combustion Instability Investigation (Ref 12) 
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IV.C, Application of the ROCCID Methodology, (cont) 


(1) Largest negative and positive growth coefficients for the most stable and 
unstable modes without acoustic damping. 

(2) Highest sensitivity of combustion stability to chamber pressure and mixture 
ratio. 

(3) Greatest design change from existing data base. 

(4) Capable of stabilization at nominal operating conditions using an acoustic 
cavity. 

(5) Highest confidence in modelling dropsize and mixing efficiency. 

(6) Ease of incorporating design features into existing hardware. 

Selection of the Injector 

Three injector designs, a fine triplet, a coarse triplet, and a like-on-like doublet, were 
selected to be evaluated by ROCCID for detail design and testing. Initially, the point design 
option ROCCID was used to size the injectors for the nominal operating conditions. Thereafter, 
the analysis portion of ROCCID was used to perform trade studies on the initial designs. The 
experience gained from the experimental results of Ref. 13 was also applied to the analysis. Also 
from past experience, certain model combinations in ROCCID were used because they produced 
good results. HLFI was used to analyze the chamber response. Smith-Reardon and N/x correla- 
tions were used to analyze the burning response. INJ was used to analyze the injector response. 

The injector orifice sizes evaluated were 0.090 in. for the fine triplet, 0.159 in. for the 
coarse triplet, and 0.100 in. LOX and 0.065 in. fuel for the like-on-like doublet. The predicted 
performance efficiencies for all three injectors were greater than 97 percent so performance did 
not become a critical concern in the selection process. 

An extensive number of stability predictions were generated and only a portion will 
be discussed in the following section. References 14 and 15 should be examined for more 
details. 


The chug stability of the three injectors was analyzed. The fine triplet and the 
doublet had low chug thresholds at about 300 psi. The coarse triplet had a more undesirable 


RPT/P0213.105-IV 


60 


9/30/92 


IV.C, Application of the ROCCID Methodology, (cont) 

chug threshold at about 600 psi. Considering that the nominal operating pressure is 1250 psi, 
this higher chug pressure limit for the coarse triplet would not allow for much leeway in 
throttling the engine to find stable and unstable high frequency stability regions. 

The high frequency stability characteristics of the three injectors were also analyzed. 
The coarse triplet was found to be the most stable injector design (?i less than -37 in all modes) 
without acoustic damping devices. From a flight hardware designer's point of view, this would 
be an excellent design. However, since our validation criteria require both stable and unstable 
operation for verification, the coarse triplet injector design was not sufficient for ROCCID 
validation. 


The doublet was found to be stable except for the first tangential mode. The doublet 
was very promising because it had a region of operation where it could be driven stable and 
unstable by changing the mixture ratio as shown in Figure 20. The doublet could be stabilized 
using acoustic damping devices that would allow validation of the models for damping devices. 
The major drawback of the doublet was that it was very similar to the design used in Ref. 13. 
Since this was a validation of the predictive capability of the models in ROCCID, using an 
injector so close to the anchored hardware made this injector a secondary choice. 

The fine triplet design was predicted to be unstable in several different modes with 
varying operating conditions. Figures 21 to 24 show ROCCID predictions of stability changes 
over the operating map. The first longitudinal (1L) mode was stable at all but lower chamber 
pressures. The first tangential (IT), second tangential (2T), and first radial (1R) modes were 
predicted unstable over wide changes in operating conditions. An evaluation of an orifice 
diameter increase from 0.09 to 0.10 in. indicated that the stability margin did not improve 
sufficiently to provide any significant advantages. Damping devices were found to be very 
effective in damping the unstable modes. In addition, this design was very different from the 
hardware described in Reference 13 in terms of number of elements and orifice sizes. 

Table 13 provides a rating of the candidate injector designs. It is obvious that the 
fine triplet injector was the best choice to validate the models contained within ROCCID. The 
fine triplet injector stability was predicted to be sensitive to chamber pressure and mixture ratio 
variations that would allow for testing the sensitivity and predictive capabilities of the models. 
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Figure 20 


Figure 21 


Figure 22, 
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Figure 25. Variation of Growth Coefficient X With N and x 
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TABLE 13 


INJECTORS COMPARED AGAINST RATING CRITERIA 


Rating Criteria 

Injector Configuration 

Fine O-F-O 
Triplet 

Coarse O-F-O 
Triplet 

Like- on- like; 
One Doublet 

Largest negative and positive 
growth coefficients for least 
damped mode without acoustic 
damping 

Host unstable 

Most stable 

Unstable 

Best sensitivity of combustion 
stability to P c and 0/F 

Significant variation 
with 21 and IT; stable 
operation with 0/F and 

Much larger 
variation in P c 
and 0/F required 
to drive unstable 

Good sensitivity with 
O/F for IT mode 

Greatest change from existing 
data base 

N„ 39 to 105 
D of 0.125 to 0.090 


N„ 39 to 48 
D 0 , 0.125 to 0.159 

N„ 105 to 123 
D„, 0.059 to 0.065 

Capable of stabilization at 
nominal operating conditions 
using an acoustic cavity 

Three types of 
cavities evaluated 


Stable without an 
acoustic cavity 

IT monotune 
cavity judged to 
be sufficient 

Best confidence in modeling 
drop size and mixing 
efficiency 

r m reasonable; 
good 


r m reasonable; 
E m good 

r m reasonable; 
E m reasonable 

Ease of incorporating design 
features into existing hard- 
ware 

New faceplate only 


New faceplate and 
injector core 

New faceplate and 
injector core 

Preferred rating 

6 of 6 

3 of 6 

2 of 6 


Indicates preferred 
characteristic 


IV.C, Application of the ROCCID Methodology, (cont) 

The fine triplet injector was sufficiently different from previously tested injectors to allow a true 
" a priori" prediction. As a result, this injector was selected as the validation test hardware. 

Growth Coefficient Uncertainty 

The sensitivity of the stability predictions (in terms of the growth coefficient) to the 
pressure interaction index, N, and sensitive time lag, z, were evaluated over a range of N and z 
values. The ranges of N and z values were selected based on estimates from three sources. 

(1) The empirical variations in these parameters shown in Refs. 16 and 17. 

(2) The calculated variation resulting from different droplet combustion mixture 
ratios and gas temperatures using the CRP model (Ref. 9). 

(3) Variations in t based on an estimated 1 percent uncertainty in energy release 
efficiency and relating that uncertainty to a change in vaporization time lag. 

The three techniques yielded similar results. Based on the Ref. 16 data the uncer- 
tainly in N was estimated to be ±26% (0.90±0.20) and the uncertainty for z was +43%, -31% (0.7 
x 10 -4 sec to 1.79 x 10" 4 sec) (see Figures 25 (c) and (d) of Ref. 16). 

The CRP model was used to evaluate assumed fuel droplet surface mixture ratios of 
0.7 to 2.4 and gas temperatures of 2500 to 6500°R at the droplet surface. These analyses resulted 
in N values of 0.45 to 1.0 and z values of 0.94 x 10' 4 sec to 2.5 x 10" 4 sec, similar to the ranges 

shown from Ref. 16 but slightly larger. 

Additionally, a typical 1% performance uncertainty (98 ±1%) resulting from a fuel 
droplet size uncertainty would yield a z uncertainty of approximately 20%, again in the general 
range observed from the Smith- Reardon data. 

The above variations in N and z from Smith-Reardon at the nominal chamber 
pressure and mixture ratio (Pc = 1250 psi, O/F = 2.8) were evaluated using ROCCID. The 
results, as shown in Figure 25 indicate large variations in predicted growth coefficient, up to 
approximately ±600 sec' 1 for the 2T mode if both N and z variations are considered. However, 
implicit within ROCCID is a variation in z due to engine operating conditions. Therefore, the 
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IV. C, Application of the ROCCID Methodology, (cont) 


data scatter leading to uncertainties in % from Refs. 16 and 17 may be due to differences in the 
operating condition of the engines being evaluated. If only the uncertainty in N is considered, ^ 

the growth coefficient in variations are reduced to ±25 sec_1 for the 1L mode, ±350 for the 
IT mode, and ±500 for the 2T mode. Based on the statistical scatter in the existing empir- , 

ical correlations of N and the present state of the models, these large bands of uncertainty asso- 
ciated with the predictions are expected. 

It should be noted that there is not an equal expectation of N values within the pub- 
lished uncertainty band of ±25%, which probably represents a ±3a standard deviation. If this is 
the case, then a statistical analysis of experimental data going into ROCCID could be performed 
to estimate the confidence level of the ROCCID predictions for calculated values of growth 
coefficient. A more detailed analytical and statistical evaluation of the data is recommended to 
allow a better understanding of the reasons for and consequences of the uncertainties in N and t. 

D. ROCCID VALIDATION TEST PLAN DEVELOPMENT 

A plan for generation of a test data base for validation of the ROCCID Methodology 
was developed that consisted of nine basic test series identified as test series A-I in Table 14. 

These test series were developed to verify the thrust chamber start and shutdown sequence (series 
A), stable operation and performance with an acoustic resonator (series B and D), the chug 
stability predictions (series C), and unstable operation without acoustic damping (series E and 
G). Test series to verify stable operation without acoustic damping (series F) and unstable 
operation with acoustic resonators (series H) were also included to assess ROCCID’s ability to 
forecast major combustion stability changes through operation at extreme off-nominal operating 
conditions. If tests with an acoustic resonator show much larger regions of unstable operation 
than predicted by ROCCID, a test series without an acoustic resonator (Series I) would be con- 
ducted to evaluate the instability frequency and amplitude suppressions (if any) resulting from 
the cavity test series to aid in cavity redesign. 

+ 

A total of twenty six tests were planned, but the total number of tests and specific 
sequence and test series were subject to change pending program cost and schedule constraints * 

and potential test results. A specific test matrix and test sequence for this testing is included in 
Table 15. The Table 15 information provides a test sequence overview. More details of the test 
logic and expected results are contained in the discussions of Block 1, Block 2, and Block 3 
testing in subsequent paragraphs. 
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TABLE 14. ROCCID VALIDATION TEST SERIES 
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TABLE 15 . ROCCID VALIDATION TEST MATRIX 
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(1) See Table VI for Test Series Description 

* Identifies Tests of Lower Priority That May Be Eliminated in the Event of Unforeseen Problems That Require Additional Tests 
To Achieve Adequate Resutls 


IV.D, ROCCID Validation Test Plan Development, (cont) 

The testing was planned in three blocks, as shown on the overall test logic diagram. 
Figure 26. Block 1 tests used a 1/4 wave tube acoustic cavity and the tests were structured to 
meet objectives of test series A (start, shutdown, and flow rate balance verification), B 
(verification of stable operation at the nominal operating point, MR = 2.8, Pc = 1250 psia), C 
(verification of chug and L mode predictions), and D (verification of stable operation at off 
nominal conditions). Block 2 tests were conducted without chamber acoustic damping devices 
and were structured to meet the objectives of test series E (verification of unstable operation w/o 
damping devices at the nominal design point), F (verify stable operation w/o acoustic damping 
devices at off nominal conditions), and G (verify unstable operations w/o acoustic damping 
devices at off nominal conditions). If the Block 1 tests did not result in a significant region of 
stable operation a retuned or re -configured acoustic cavity would be designed, fabricated, and 
tested based on the Block 1 cavity temperature and stability data. In this case, the Block 2B 
testing (without a cavity) would be conducted to confirm undamped chamber modes at the 
nominal operating condition and over a range of flow rates. These tests would be used to assess 
the undamped configuration ROCCID predictions and the mode suppression observed in the 
Block 1 tests would further verify the Block 1 cavity sound speed. The tests would be conducted 
while the retuned cavity is being prepared for testing. 

Block 3 testing, with an acoustic cavity, was used to confirm nominal operating point 
stability and verify predicted off nominal operating point unstable operation if, based on Block 1 
tests, the original cavity tune was adequate. If the original cavity tune was not adequate to 
provide a substantial region of stable operation based on Block 1 tests, Block 3B tests would be 
used to evaluate a better cavity tune or configuration, using the ROCCID results anchored to the 
Block 1 data. 

Tests identified with an asterisk in Table 15 are designated as lower priority than the 
remaining. In the event that additional unplanned tests were required to achieve any specific test 
objective, these tests would be considered for deletion from the test program in order to maintain 
control of the overall program cost. 

As will be shown, each of the three blocks of tests contain 7 to 10 tests. Based on 
previous experience, this would allow one chamber liner to be used for each block of tests. 
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IV.D, ROCCID Validation Test Plan Development, (cont) 


Block 1 Tests 


More detailed test logic along with specific success criteria is shown on Figure 27 
for the Block 1 tests. One should note that the generation of well defined test data, i.e., accurate 
flow rate, temperature, pressure and thrust data (and for the unstable tests accurate combustion 
frequency determination) were of paramount importance. These data were used not only for 
comparison with ROCCID model predictions but also provided the basis for updating or re- 
anchoring the model, should that be necessary. 

The specific Block 1 test points are shown tabulated on Table 16 and shown graph- 
ically on the test facility operating map on Figure 28. The predicted regions of stable and 
unstable operation shown on Figure 28 were based on the ROCCID growth coefficient 
predictions. 


The test series A logic is shown on Figure 29. Test series A verified the start and 
shutdown sequencing, the flow system resistances and the instrumentation (except CSMs) before 
proceeding onto the performance and stability tests. Two tests were allocated for this test series. 

The logic for test series B, verification of nominal operating point stability, is shown 
on Figure 30. This was the first test of sufficient duration to verify cavity temperature estimates, 
provide meaningful performance data and the first test to assess bomb induced chamber pressure 
spikes. Four important “gates” are shown on the test logic diagram. The first gate is related to 
the success criteria of Figure 27 (Is the test data valid and was intended operating point 
achieved?) If the answer was yes the next gate would be evaluated. If the answer was no, the 
test data and prediction would be evaluated and a decision to retest or proceed to next gate made. 

If the test was statistically stable and the bomb overpressure is less than 20% of 
chamber pressure, a retest would be conducted. If the test was spontaneously unstable, the mode 
was defined and a retest made at a more stable operating point defined from the growth coef- 
ficient predictions based on the ROCCID analysis. 

The test series C logic, shown on Figure 31, was used to ensure that valid data was 
collected and evaluated for verification of chug and longitudinal mode instabilities. The 
ROCCID prediction of longitudinal mode instability was the result, to a large extent, of 
predictions of injection coupling. This can be seen from the response plots for MR = 2.8 and 
chamber pressures of 625 and 800 psia shown on Figure 32. Injection coupling with impinging 
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Figure 27. Block 1 Test Logic and Success Criteria 
(Chamber With 1/4 Wave Damping) 


TABLE 16. BLOCK 1 TESTS ARE CONDUCTED WITH AN ACOUSTIC CAVITY 
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Figure 28. Block 1 Tests Will Define Stable Operating Region With an Acoustic Cavity and 
Evaluate Chug and L Mode Instabilities 


OBJECTIVE 

Verify Start System, Flow Systems and Test Instrumentation 
Before Proceeding to Performance and Stability Tests 
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Figure 29. Block 1 - Test Series A Logic 
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Verify Stable Operation at Nominal MR-Pc With 1/4 Wave 
Cavity/ Define and Verify Stable Operating Point 
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Figure 30. Block 1 - Test Series B Logic 


OBJECTIVE 

Verify Chug and L Mode Stability Predictions 
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Figure 31. Block 1 - Test Series C Logic 
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Figure 32. Stability Predictions at Low Chamber Pressure 
Show Large Influence of Injection Coupling 
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IV.D, ROCCED Validation Test Plan Development, (cont) 

stream injectors is not well documented, and thus the comparison of the model prediction with 
the test data provided valuable ROCCED calibration data. 

The final Block 1 test series, test series D, was structured to define the operating 
region of stable operation with an acoustic cavity. Test series D logic is shown on Figure 33. 

The series D tests were to be conducted in the order shown on Figure 28 if, as predicted, all the 
tests were stable. If an unstable test occurred, the mode of the instability and the appropriateness 
of the cavity tune (based on measured cavity gas temperature) would be defined and the growth 
coefficient data used to determine the most stable operating point for that mode. If a re-test at 
this operating point was stable, the next untested test point from Table 16 would be tested. If the 
re-test is unstable in the same mode as the previous test a second re-test at a maximum or mini- 
mum throttle condition would be conducted to determine the extent of the unstable region for that 
mode. If the re-test was unstable in a different mode from the original test, a new test point 
would be selected for the new mode from the growth coefficient data. 

Block 2 Tests 

As shown on the logic diagram of Figure 26, the Block 2 testing was conducted 
without an acoustic cavity. Block 2A tests were conducted to verify unstable operation without 
an acoustic damping device at nominal and off nominal operating conditions (test series E and G 
respectively) and verify stable operation at extreme off nominal conditions, test series F. The 
Block 2A test series logic is shown on Figure 34. Specific test conditions for these tests are 
shown on Table 17 and further illustrated on the test facility operating map on Figure 37. 

As shown on the Figure 35 logic diagram and the specific test listing of 
Table 17, the first test of this series (Test #11) was at nominal mixture ratio and chamber 
pressure. This test was predicted to be unstable in the IT or 2T mode. If this was verified by the 
test data the next test was at a predicted stable condition as shown on Table 17. If Test 1 1 is 
stable and the success criteria of Figure 27 is met, one re-test would be conducted. If this re-test 
was also stable, the ROCCED predictions were used to define the maximum unstable operating 
point and this condition was tested next. If subsequent tests were also stable the process was 
repeated until all the operating points with the maximum predicted growth coefficient for each 
mode (most unstable predicted condition) were surveyed. This resulted in a survey of a wide 
range of test conditions that was used for ROCCID model anchoring. If tests that were predicted 
to be stable were found to be unstable, another predicted stable point was tested. If a predicted 
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Figure 33. Block 1 - Test Series D Logic 
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Figure 34. Block 2A Test Series Logic 
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Figure 35. Block 2A Tests Will Define Stable and Unstable Operating Regions Without an 
Acoustic Cavity 


IV.D, ROCCID Validation Test Plan Development, (cont) 


stable point was in fact stable, the test was repeated to verify the stable result.Tests 17-19 were 
allocated for these tests. 

If, as a result of Block 1 testing, the stable operating region with an acoustic cavity 
was much smaller than predicted, a retuned (different cavity depths) or reconfigured (different 
1T/2T cavity numbers or % area) acoustic cavity would be made based on the Block 1 test data. 
In this event, the Block 2 non-cavity testing would be used to determine the amount of frequency 
and amplitude suppression that was achieved by the cavities by comparing the Block 1 data to 
data collected during the Block 2B test series. This test series is shown on Table 18. 

Block 3 Tests 


As shown on the logic diagram of Figure 26, the Block 3 tests would have 
reconfirmed the stable operating region and defined the unstable operating region with the 
original cavity tune from the Block 1 tests (Block 3A); or as actually happened if a new cavity 
configuration was required to provide stable operation, and the Block 3B tests were used to 
evaluate the new cavity design. The Block 3 test logic is shown on Figure 36 and the specific 
test points are shown on Table 19 and Figure 37. 

Instrumentation Requirements and Accuracy 

The objective of thrust chamber testing was to generate a data base for valida- 
tion of performance and combustion stability predictions made using the ROCCID methodology. 
Therefore, the test hardware was heavily instrumented to provide the required data base. 
Instrumentation was also included to control and monitor the facility operation and thrust 
chamber transients and provide for test kills in the event that red lines were exceeded for 
specified critical parameters. 

The data acquisition requirements for determination of combustion stability and per- 
formance during this testing are provided in Tables 20 and 21, respectively. Parameter ranges 
and accuracy requirements are specified on these tables. In some cases a single measurement 
was used to derive more than one parameter resulting in different measurement accuracies. 

A complete list of instrumentation, range and recording device requirements was 
provided in the test plan (Ref. 1 8). The location of facility instrumentation is shown in the test 
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TABLE 18. Block 2B Tests Were Conducted Without an Acoustic Cavity to Provide 

Data for Acoustic Cavity Re-Design 
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TABLE 1 9. BLOCK 3 TESTS WERE MADE TO CONFIRM STABLE OPERATION WITH AN 
ACOUSTIC CAVITY 
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Figure 37. Block 3A Tests Were Planned to Re-Confirm Stable Operating Regions With an Acoustic 
Cavity and Define Regions of Unstable Operation With the Cavity 


TABLE 20. COMBUSTION STABILITY DATA ACQUISITION REQUIREMENTS 


1) Dynamic Chamber Pressure 

• Measurement 

Frequency: 100 to 10K hz 

Amplitude: ±100% Pc 

Damp & Growth Rates: psia/milliseconds 

• Location 

Minimum: 3 clock positions near 
the injector face and 3 axial 
locations at the same clock position 

2) Dynamic Manifold Pressure 

• Measurement 

Frequency: 100 to 5K hz 

Amplitude: ±33% Pmanifold 

• Location 

As close to orifice inlet as feasible 

3) Operating Conditions 

• Static Chamber Pressure: ±2% accuracy 

• Mixture Ratio: ±3% accuracy 

• Propellant Temperature: ±10°F 

• Static Manifold Pressure: ±2% 

4) Combustion Perturbation (Dynamic Rating) 

• Instantaneous 

• Minimum 25% overpressure 

• Close to the injector face 

5) Acoustic Cavities (if present) 

• Gas Temperature: ±100°F 

6) Accelerometer 

• Attached to the thrust chamber flange 


89 


Defined Consistent 
with CPIA 148, 170 
and 247 
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TABLE 21. THRUST CHAMBER PERFORMANCE DATA ACQUISITION 

REQUIREMENTS 


1) Vacuum Specific Impulse ±1% 


Thrust Level: 
LOX Flowrate: 
RP-1 Flowrate: 


I 


Sufficient to provide overall 
specific impulse 
accuracy of ±1% 


2 ) 


Static Chamber Pressure: ±1% 

Mixture Ratio: ±2% 

Tolerance Requirement is usually 
Propellant Temperature: 

Expansion Area Ratio: 

Ambient Pressure: 

Nozzle Exit Diameter: 

Inert Purge Flow: 

Characteristic Velocity ±2% Accuracy 
Chamber Pressure ±1% accuracy 


set by flowrate accuracy requirements 


±5°F 

± 2 % 

±0.1 psia 
±.05 in. 

±.05% of Wj Reactive 


LOX Flowrate: 
RP-1 Flowrate: 


£ 


Accuracy requirements set by 
specific impulse treatment 


Nozzle Throat Diameter: ±.02 in. 
Mixture Ratio: ±2% 


Propellant Temperature: ±5°F 

Inert Purge Flow: ±.05% of Wj Reactive 

3) Manifold Pressure: ±1% 

4) Static Chamber Pressure Profile 

Measurements sufficient to determine the total pressure loss from the injector face to the 
start of nozzle convergence. An array of static pressure transducers and/or AP gauges 
with in-place calibration to measure AP (station-to-station) of 5 to 25 psia. 
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IV.D, ROCCID Validation Test Plan Development, (cont) 

facility schematic was included in the test plan. Location of the thrust chamber instrumentation 
was noted on the hardware drawings. 

E. VALIDATION HARDWARE DESIGN 

The hardware design for the test validation of the ROCCID methodology consisted 
of a thrust chamber that operated at moderate chamber pressure (1,250 psia) and thrust (48,000 
Lbf) using LOX/RP- 1 propellants. The thrust chamber was designed for short duration operation 
at sea level conditions adequate to experimentally determine its performance and combustion 
stability over a wide range of chamber pressure and mixture ratio. The drawing tree shown in 
Figure 38, identifies the components and their Aerojet drawing numbers that made up the engine 
assembly for ROCCID validation testing. A complete collection of these drawings is provided in 
Appendix D of Volume II of this final report and in Reference 19. 

The basic hardware components were originally designed, fabricated and tested 
during the Oxygen/Hydrocarbon Injector Characterization Program (Contract F0461 1-85-C- 
0100) conducted over the time period from September 1985 to January 1990. The critical design 
features of this hardware, such as the injector element faceplate, acoustic cavities and the com- 
bustion chamber internal (gas-side) contour including throat size, contraction ratio and chamber 
length were redesigned using the ROCCID methodology. These changes were incorporated into 
a revised design disclosure (drawing package) and the hardware components were appropriately 
modified and refurbished for this test program. This section includes a description of this thrust 
chamber and its components. 

1. Engine Assembly 

The engine assembly for use in this test program is shown in Figure 39. Major 
components of this thrust chamber assembly include: (1) an injector assembly consisting of an 
outer flange and RP-1 manifold and an inner core with a liquid oxygen inlet and the OFO triplet 
injector faceplate, (2) an ablative (silica phenolic) lined combustion chamber with a high strength 
steel (4340) shell designed for operation at chamber pressure as high as 1,750 psia, and (3) ancil- 
lary components including a proof and leak check plate, seals, bolts, nuts and washers, bomb 
adapters, and acoustic cavity inserts. 

This thrust chamber was designed to interface with the propellant feedsystem 
and thrust takeout existing in the Aerojet Propulsion Division's E-4 test stand. The fully 
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Figure 38. The Drawing Tree for the ROCCID Test Validation Thrust Chamber 
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Figure 39. The Engine Assembly for ROCCID Test Validation 


IV.E, Validation Hardware Design, (corn) 


assembled thrust chamber weighs an estimated 1400 lbm. A photograph of a fully assembled 
thrust chamber on the E-4 test facility with instrumentation in place is provided in Figure 40. 

The hardware design included all ports and fittings for the required instrumen- 
tation. Ports for 14 static chamber pressure transducers, 6 dynamic chamber pressure trans- 
ducers, static and dynamic pressure transducers in each propellant manifold and 4 thermocouples 
to measure acoustic cavity gas temperatures were included. In addition, two ports were provided 
in the combustion chamber for the bomb adaptors that contain RDX explosives to induce 
chamber overpressure spikes for measurement of dynamic combustion stability characteristics. 

2. Injector Assembly 

A drawing of the injector assembly (P/N 1206427) is shown in Figure 41. The 
assembly consisted of a flange assembly (P/N 1200773-19) shown in Figure 42 and an injector 
core assembly (P/N 1206423) shown in Figure 43. The flange assembly contained the fuel inlet 
and the thrust takeout and chamber mounting surfaces. The injector core assembly contained the 
oxygen inlet, the fuel and oxygen downcomers, distribution channels and the injector faceplate. 
The injector core assembly was designed to bolt to the injector flange assembly. An assemblage 
of photographs showing the injector assembly and its components fabricated for previous 
LOX/Hydrocarbon testing is included in Figure 44. 

The injector flange assembly was made from 304 stainless steel and weighs an 
estimated 430 lbm. This assembly was attached to the test facility using sixteen (16) 1/2 inch 
bolts. A 2 inch diameter Greyloc (GR 20, schedule 160) fitting was included for connection to 
the test facility fuel delivery system. 

The injector core assembly, shown in Figure 43, consisted of an injector core 
made from 304 stainless steel, a schedule 160 oxidizer inlet tube with a 3 inch diameter Greyloc 
(GR 25, schedule 160) fitting, a zirconium copper faceplate and a stainless steel igniter line 
assembly. The injector assembly inserted into the flange assembly and was attached using 16 - 
5/8 inch bolts. Two teflon O-ring seals provide a seal between the flange and injector core 
assemblies. 


The injector core included a zirconium copper faceplate that was brazed to the 
304 stainless steel body. The faceplate contained the fine thrust-per-element OFO triplet injector 
pattern that was developed using the ROCCED methodology. The 105 OFO triplet elements 
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Figure 40. The ROCCID Test Validation Engine Assembly Installed on Test Stand E-4 



Figure 41. The Injector Assembly Used an Existing Fuel Manifold 
and Injector Core and a New Injector Faceplate 
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Figure 42. The Injector Flange and Outer Fuel Manifold was Available 
From a Previous Phillips Laboratory Program 
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Figure 43. An Existing Injector Core Was Used 
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Figure 44. The LOX/Hydrocarbon Injector Assembly is Test Proven 


IV.E, Validation Hardware Design, (corn) 


were contained in three circular rows around a centrally located TEAL/TEB ignition port. Each 
OFO injector element consisted of two .090 in. diameter orifices that inject liquid oxygen on a 
central RP-1 stream formed from a central axially directed center .090 in. diameter orifice. A 
drawing of the injector faceplate showing element design details is provided in Figure 45. The 
.090 in. orifices were the smallest fluid flow passages in either the liquid oxygen or the RP-1 
propellant circuits. 

The design details of the O-F-O triplet element are provided in Figure 46. The 
two oxygen orifices were canted at 35° with respect to the axially directed central fuel orifice. 
The three orifices were aligned along a radius of the injector and spaced so that their centerlines 
impinge at a single point 0.288 in. off the face of the injector. The faceplate was 0.312 in. thick 
resulting in an orifice L/D of 3.5 for the fuel orifice and 4.3 for the oxygen orifices. The inlets to 
all orifices were countersunk with a 0.1 15 in. diameter that provided an orifice to inlet area ratio 
of 0.61. The backside of the copper faceplate was electropolished to round all surfaces on the 
orifice inlet. Cold flow tests (Ref. 13) have shown that an orifice Cd of .91 to .93 is achieved 
with this type of orifice geometry. 

The TEAL/TEB igniter circuit contained in the injector core consisted of a 
flared 1/4 in. OD CRES tube that attached to the test facility using a nut attached to the tubing. 
The igniter runs through the center of the injector core and the TEAL/TEB was injected through 
a CRES igniter plug that contained 3-. 124 in. ports. 

3. Thrust Chamber Assembly 

A drawing of the combustion chamber assembly is shown in Figure 47. The 
assembly consisted of a high strength 4340 steel structural shell and a single piece replaceable 
silica phenolic ablative liner. The weight of the thrust chamber assembly was estimated as 
approximately 600 lbm without the closure plate and instrumentation. 

The chamber assembly was attached to the injector assembly with 20 1-inch 
diameter bolts. The retainer plate was attached to the steel chamber housing using 16 1-inch 
diameter bolts. O-ring seals were used at both the injector/chamber and chamber/retainer plate 
flanges to provide a hot gas seal. 

The hot gas combustion chamber internal profile was formed using a silica 
liner, shown in Figure 48 that mated with the steel shell. The liners were inserted into the steel 
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Into a Zirconium Copper Faceplate 




<o 


102 


Figure 46. The OFO Triplet Element was Designed to ROCCID Specification 
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Figure 47. A New Thrust Chamber Design was Prepared to Meet 
ROCCID Validation Requirements 
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Figure 48. A New Silica Phenolic Liner Design was Created to Match ROCCID Design Criteria 


IV.E, Validation Hardware Design, (cont) 


shell and bonded in place using RTV. A retaining ring (P/N 1200976), was also used at the aft 
end of the combustion chamber to locate and support the liner within the chamber shell. 

The replaceable silica phenolic liner was designed with a throat diameter of 
5.43 in. and a cylindrical combustion chamber with a diameter of 7.68 in. This resulted in a 
contraction area ratio of 2: 1 (area of the nozzle inlet to the nozzle throat). A conical expansion 
nozzle having a 4.5° half angle was used to expand the combustion gases to a 1.6:1 expansion 
area ratio. The liner was designed to provide a chamber length (injector face to throat plane axial 
distance) of 14 in. 

Previous test experience has shown that for short duration (e.g., 0.3 sec steady 
state) stability and performance tests that the silica phenolic material is superior in terms of 
durability to either the bulk graphite or carbon phenolic liners. The carbon based liners are sus- 
ceptible to oxidation, particularly with the OFO injector pattern without fuel film cooling. Based 
on this test experience, the liners were expected to last for 6 to 10 tests before requiring 
replacement. 


Acoustic damping in the form of an acoustic resonator was supplied using a 
replaceable OFHC copper ring that was installed at the head end of the combustion chamber. A 
bituned quarter wave acoustic cavity design using the ROCCID methodology is shown in 
Figure 49. The cavity tune was adjusted by changing the depth of the slots within the copper 
ring. A solid copper ring (without slots) was used when no acoustic damping was desired. Refer 
to Figure 39, which shows the installation of the cavity block (item (D and ®) with and without 
acoustic cavities. 


Bombs, which were used to perturb the combustion process to measure 
damping characteristics, were held in the combustion chamber using an adapter (P/N 1201080) 
shown in Figure 50. The adaptor body was fabricated from 1040 carbon steel and was threaded 
into the steel chamber shell in two locations. The adapter contained a carbon phenolic insert and 
teflon bomb holder that extended through the chamber liner to a position slightly recessed from 
the gas side liner surface. This recess was filled with RTV 60 to provide a measure of thermal 
protection to the bomb during hot-fire to preclude thermal detonation. This insert was 
replaceable and was designed to shield the chamber liner from the effects of the bomb. When no 
bomb was used, this insert was replaced with a carbon phenolic plug which filled the entire bomb 
cavity within the chamber liner. 
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Configuration With A Bomb Configuration Without a Bomb 








107 


Figure 50. Bomb Adapter Design With a Thermal Protection Enhancement 


IV, Technical Narative, (cont) 


The chamber contained provision for a multitude of instrumentation as shown 
on Figure 51. Included were 6 high frequency pressure transducer ports, 14 static pressure ports, 
two bomb ports, and 4 gas temperature pons. These ports were supplied in the steel shell and 
drilled into the replaceable silica phenolic liners after assembly. 

F. HARDWARE FABRICATION 

Component fabrication in support of the planned ROCCID hot fire test validation 
program was successfully completed over the period June 1990 through March 1991 (Reference 
20). A hardware inventory necessary to support the approved test plan (Ref. 18) was obtained 
through rent free use of thrust chamber components from the AF 0461 1-85-C-0 100 contract and 
from purchased parts and hardware fabrication subcontracts. 

A complete list of components required for the ROCCID Validation Thrust Chamber 
Assemblies (P/N 1206430-9, -19, and -29) is provided in Table 22. Spare parts were included to 
enhance the probability of successfully completing the test program in a timely and cost efficient 
manner. For example, two identical injector cores and two chamber shells were included. The 
backup injector core, shown in Figures 52 and 53 was included as a risk reduction item because it 
was the most vulnerable component from the standpoint of survival of the combustion 
instabilities that were expected during validation testing. Two chamber shells were included so 
that replaceable liner installation could be performed on one chamber while testing was 
conducted with the other. Four silica phenolic ablative liners were provided for the planned 26 
tests. Three acoustic resonator block assemblies were included: one machined to provide the 
specified 1/4 wave tube bituned cavity (P/N 1206426) and two assemblies that totally blocked 
the acoustic cavity region. Unfortunately, a design drafting error in the acoustic cavity blocks 
(PN 1206426 or 1206431) went undetected through both the design and test phases resulting in 
acoustic cavity configurations (both "open" and "blocked") that were far from the desired design 
as noted in Figures 54 to 56. 

The acoustic cavity was formed from the assembly of five hardware components (the 
injector core, flange, chamber, liner and the acoustic block), and therefore the actual cavity and 
its critical dimensions that affected its "tune" were not shown on the component drawings. 
Therefore, inspections of individual drawings were not sufficient to reveal this design flaw. The 
effect of the cavity design errors on the operation of the combustion chamber and the pretest 
combustion stability characterization is discussed in Section VI. H of this report. 
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Figure 51. Provision for Instrumentation Required for Validation 
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Figure 52. The Injector Core With Its Faceplate Installed for Brazing 
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Figure 54. Bitune-Cavity Configuration Comparison of Actual and Intended Acoustic Cavities 
For ROCCID Testing 
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Figure 56. Monotune Cavity Configuration 


IV. F, Hardware Fabrication, (cont) 


Three types of verification tests were conducted to determine the adequacy of the 
hardware for hot fire testing. 

First, the integrity of the copper faceplate braze joint to the stainless steel injector 
core was verified through a moderate pressure leak check using helium. During this test, the 
injector faceplate was sealed using a soft rubber mat and metal plate held over the face using C- 
clamps. The fuel circuit was filled with water and the oxygen circuit was pressurized to approx- 
imately 50 psig and held for approximately 5 minutes. No observation of helium bubbles in the 
water filled fuel circuit verified that the braze joints of both injector cores were leak free. 

Second, all pressure containing hardware components were subjected to a proof test 
using filtered water at 3800 ±100 psig and held for 5 minutes. Components that were previously 
verified through proof testing during the F0461 1-85-C-0100 contract did not require further 
proof verification. All proof tests were successfully conducted in the Aerojet E-Zone test facility 
and were witnessed by the test engineer. 

Last, cold flow tests of both injector cores were conducted to verify flow resistance 
and the element impingement. These tests were conducted in the A-Area facility and the 
Production Hydrolab facility over a pressure and flowrange that precluded orifice cavitation 
(below approximately 50 psig). A photograph of the refurbished injector core at very low AP (3- 
5 psig) is shown in Figure 57. 

The impingement process of both injector cores was visually inspected by the project 
engineer and the lead combustion analyst. In some cases, poor element impingement was noted, 
which was subsequently attributed to blocked, or partially blocked orifices resulting from inad- 
vertent contamination from the manufacturing process. In all cases, satisfactory impingement 
was achieved after backflushing to remove the contamination. 

Non-cavitating flow resistance (R) or admittance (Kw = w /(AP Spg) 1/2 ) values were 
also measured for both injector cores as summarized in Table 23. 
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Figure 57. Water Cold Flow Tests Were Conducted to Verify 
Element Impingement and Flow Resistance 


IV. F, Hardware Fabrication, (cont) 


TABLE 23. ROCCID INJECTOR COLD FLOW ADMITTANCE 
VALUES (lbm/sec/(psid) 1/2 ) 


Injector Core 

Fuel Circuit 

Oxygen Circuit 

-0100 Refurbed New ROCCID 

2.85 - 3.0 

5.8 -5.9 

New Injector Core 

2.95 

5.7 - 5.8 


These Kw values were within the expected range thus indicating that the desired high 
orifice Cd (> .9) was achieved. 

Minor design refinements were made during the fabrication phase to facilitate the 
manufacturing ease and correct errors and omissions on the drawings. These changes have been 
incorporated into the design disclosure package using Aerojet’s Advanced Drawing Change 
Notice (ADCN). The ADCN’s prepared in support of this program are provided in Appendix E, 
of Volume II of this final report. 

In addition to these design changes some features of the as-fabricated parts did not 
meet drawing requirements. In each case, these discrepancies were evaluated by the Aerojet 
Engineering Review Board for this program to determine that they did not affect the form, fit or 
function of the part. A description of the discrepancy and its disposition are recorded on Aerojet 
Nonconformance Reports that are included in Appendix F of Volume II of this report. 

G. VALIDATION TESTING 

Hot fire tests, in accordance with the approved test plan (Ref. 1 8), were conducted on 
the Aerojet E-4 test stand in the spring of 1991. The purpose of this testing was to provide a test 
data base with hardware designed using the ROCCID methodology for assessment of the validity 
and capability of the methodology and identification of its shortcomings to support future 
improvement efforts. A total of 40 test firings were made, from which 27 tests produced results 
suitable for correlation with ROCCID combustion stability and/or performance predictions. 

The testing consisted of ”3 blocks" of tests. Block 1 tests were with a bituned 
acoustic cavity and were structured to provide test data in the most combustion stable engine 
configuration. These tests established performance values and mapped the regions of stable and 
unstable combustion. They were predicted to be the safest tests from a hardware damage stand- 
point since high amplitude combustion instabilities were not likely. This provided the most 
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IV. G, Validation Testing, (cont) 


benign environment to establish start and shutdown sequences and performance values. The 
Block 2 tests were without any acoustic cavity and provided a direct assessment of the benefit of 
the cavity used in the Block 1 tests. The block 3, and final group of tests provided test data with 
a monotuned acoustic cavity. The valid tests in these three blocks are shown in Table 24. 

An overview of the test logic for each of the Blocks of tests is shown on Figure 58. 
The path depicted by the shaded boxes was the logic path followed during the test program. The 
Block 1 tests with the bituned cavity were structured to meet the objectives of test series A (start, 
showdown and flowrate balance verification), B (verification of stable operation at the nominal 
operating point, MR=2.8, Pc=1250 psia), C (verification of chug and L mode predictions), and D 
(verification of stable operation at nominal conditions). The Block 2 tests without acoustic 
cavities were structured to meet the objectives of test series E (verification of unstable operation 
without acoustic cavities at the nominal design point), F (verify stable operation without acoustic 
cavities at off nominal conditions), and G (verify unstable operation without acoustic cavities at 
off nominal conditions). The Block 3B tests with a reconfigured (monotuned) acoustic cavity 
were conducted since the original bitune cavity did not provide as large a region of stable oper- 
ation as originally desired. This block of tests met objectives B (confirm nominal MR-Pc oper- 
ation with a new cavity tune), D (confirm off nominal stability with a new cavity), and H 
(confirm off nominal unstable points with a new cavity). 

Table 25 is a summary of the test data. It includes directly measured data, quantities 
calculated from measured data, and brief comments. Tests for checking out the engine operation 
and unsuccessful tests were not included in the table. Plots of measured system pressures for 
these tests as a function of time are included in Volume II of this report. 

The raw data was divided into two segments, performance data and liner data. The 
performance data includes the data summary period, system pressures and temperatures and 
thrust. If the test was of sufficient steady-state duration (>200 msec), the listed data is the period 
with the least variation. For the shorter tests, the data summary period typically is a single data 
point, and it is usually chosen to correspond with the time of maximum valve opening position. 

It should be noted that the flow rate values in many cases were obtained from graphical 
flowmeter records, since the digital time averaging tended to yield erroneous values because of 
transient effects in the short duration tests due to the relative location of the flowmeters and the 
thrust chamber. 
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TABLE 24. PLANNED AND ACTUAL TEST CONDITIONS 



Planned Conditions 

Actual Conditions 



Chamber 


Chamber 


Test 


Pressure 

Mixture 

Pressure 

Mixture 

No. 

Block 

psia 

Ratio 

psia 

Ratio 

4 

I 

1250 

2.80 

1177 

2.59 

7 

I 

1500 

1.50 

1441 

1.45 

9 

I 

500 

2.80 

505 

2.94 

12 

I 

300 

2.80 

250 

1.01 

13 

I 

1250 

1.20 

1210 

1.13 

14 

I 

1500 

2.80 

1358 

2.67 

15 

I 

1000 

2.00 

969 

1.93 

17 

I 

800 

7.50 

792 

6.74 

18 

I 

1250 

2.80 

1165 

2.35 

19 

I 

1250 

2.20 

1220 

2.06 

20 

I 

1250 

1.50 

1193 

1.50 

21 

II 

1250 

2.80 
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IV. G, Validation Testing, (cont) 


Liner data consists of liner number and the throat and exit diameters measured before 
and after each test. These values are typically averages of 3 measured values. 

The comment section includes pertinent comments about the test, including failure of 
key instruments, quality of performance data, and visual changes in the liner condition. The 
comments are frequently about the failure of the fuel flowmeter, that was caused by interference 
from the camera sequencer and was most prevalent at low fuel flow rates. 

The calculated data includes propellant flow rates, injector and total system flow 
admittances (Kw's) and combustion efficiencies. Propellant flow rates were calculated using the 
reported equivalent water flow rates and the propellant specific gravities, which are function of 
both temperature and pressure. The pressures at the flowmeters were assumed to be equal to the 
tank pressure less 75 psi, which is estimated to be the line loss. 

The pretest throat and exit diameters were used to calculate the throat area and the 
contraction and expansion ratios for performance calculations. The throat stagnation pressure is 
the average of the values computed using two techniques. The first technique utilizes the last 
static pressure measurement, the contraction ratio and an assumed ratio of specific heat of 1.135, 
to calculate an isentropic throat stagnation pressure. The second technique utilizes a correlation 
developed by J. Ito of Aerojet: 

Po, throat = Pc - 0.25*Pc/CR 018 

where CR is the chamber contraction ratio and Pc is the chamber static pressure. 

C* and ISP-based energy release (ERE) efficiencies were calculated as the ratio of 
the delivered to the perfect injector estimated values. In both cases, the perfect injector values 
were calculated using the TDK/BLM computer program TDK (Ref. 21) with kinetic, boundary 
layer and divergence losses from the one dimensional equilibrium (ODE) values accounted for. 

The "delivered" C*'s were computed from the throat stagnation pressures, throat 
areas, and total propellant flow rates. The throat area used was the pretest measured value. Test- 
to-test throat diameter changes were small, less than 5% from first to last test for each of the 
three chamber liners used during the testing. The delivered vacuum ISP's were calculated using 
the measured thrust, propellant flow rates, nozzle exit area and ambient pressure. 
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IV. G, Validation Testing, (cont) 


Calculations using techniques different from that described above are noted at the 
bottom of Table 25. Most common of these notes is that one or both of the flow rates were taken 
from the graphical data. In a few cases where the fuel flowmeter was interfered with by the 
camera sequencer, the mean injection Kw was used to estimate a fuel flow rate, and thus mixture 
ratio. The injection and total Kw's values are fairly consistent between tests as can be seen in 
Figures 59 and 60. The fuel and oxygen injection Kw were found to be approximately 2.94 and 
6.01, respectively. Good test-to-test correlation of the Kw’s indicated that the measurements of 
propellant flow rates, and the resulting mixture ratio were consistent. 

Figure 61 shows the PC and MR for each of the tests. The C* and ISP based effi- 
ciency values from tests that achieved sufficient steady-state duration to allow good 
measurements are plotted against mixture ratio on Figure 62. In most cases C* and ERE 
efficiencies agree within 1%, and greater agreement could be achieved by iterating on the throat 
diameter, rather than just using the pretest dimension, and by including a throat Cd term. Some 
of the scatter may also be due to lumping tests of markedly different chamber pressure into the 
same plot, although the influences of chamber pressure on the efficiencies were expected to be 
small. 


The measured cavity temperature data is plotted as a function of mixture ratio in 
Figure 63. The temperature is the average value of the temperatures measured at three different 
locations in each test. The temperature data, including values that are still in transient in several 
tests, were scattered as can be seen from the figure. As a result, these data were not sufficient to 
provide reasonable accurate speed of sound estimates inside the cavities. 

All tests were chug stable at the steady-state operating conditions. In tests 25, 28, 
and 38, chug instabilities with a small amplitude oscillation at 500 Hz occurred during the early 
portion of the transient start-up where chamber pressures were very low. These low frequency 
oscillations lasted approximately 15 to 40 msec and quickly disappeared as the chamber pres- 
sures rose to steady-state values. During the later part of the transient start-up periods, acoustic 
instabilities with higher amplitude and higher frequencies began to appear. 

A summary of the stable-combustion test results is provided in Table 26 and a sum- 
mary of the unstable-combustion test results is provided in Table 27. Both tables list test 
number, test series number (corresponding to Table 25) inferring a particular cavity configuration 
used, chamber and manifold pressures, propellant flow rates, and injected mixture 
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Figure 59. Measured Injector Kw Values Were Consistent Between Tests 
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Figure 60. Total Kw Values Were Consistent Between Tests 
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Figure 61. Test Operating Conditions 
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Figure 62. Calculated Performance Efficiency for Tests With Valid 
Steady-State Measurements 
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Figure 63. Average Acoustic Cavity Gas Temperature Measurements 
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TABLE 26 . SUMMARY OF STABLE-COMBUSTION TEST RESULTS 
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TABLE 27. SUMMARY OF UN STABLE-COMBUSTION TEST RESULTS 
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IV. G, Validation Testing, (cont) 


ratio. In addition, both tables list the high-frequency pressure transducer number from which the 
presented stability data was obtained. The propellant injection pressure drops, calculated from 
the measured manifold and chamber pressures, are also listed in the tables. In addition to the 
operating parameters, the stable-combustion table (Table 26) lists the bomb over-pressure, damp 
time, and decay (or negative growth) coefficient. These values were obtained or calculated from 
the pressure oscillograms. In addition to the operating parameters, the unstable-combustion table 
(Table 27) lists the observed resonant frequencies, growth coefficients, and amplitudes of the 
chamber pressure oscillations. These values were obtained from plots showing the evolutions of 
amplitudes and frequencies and from the power spectral analyses of the chamber pressures. 

All combustion instabilities were spontaneous. Although combustion perturbation 
bombs were used and generated over-pressures in the range of approximately 5 to 67% of the 
mean chamber pressure, no tests were driven unstable by the bombs. Power spectral density 
(PSD) analyses of chamber pressures in all unstable combustion tests showed several resonant 
frequencies existing simultaneously. The first resonant frequencies in most of the unstable com- 
bustion tests correspond to the IT acoustic frequency of the combustion chamber, which is 
approximately 3300 Hz. Other resonant frequencies appeared to be the harmonics of the IT 
mode although they may correspond to the frequencies of higher mixed mode. The resonant fre- 
quency in test number 40 observed during the early portion of steady-state operation closely 
matches the 2T acoustic resonant frequency of the chamber. 

Comparison of the chamber pressure oscillation amplitudes between the test series 
(Table 26) shows stabilizing effects of the acoustic cavities. The amplitudes and growth coef- 
ficients of the chamber pressure oscillations were generally larger in test series 2 (blocked 
cavities) than those in the other test series (with monotuned or bituned acoustic cavities). This 
indicated that the bituned and monotuned cavity configurations did provide damping to the 
system, although it was not adequate to completely suppress the instabilities. The fact that most 
of the combustion instabilities were in a first tangential mode even in test configurations with 
monotuned and bituned acoustic cavities suggested that the acoustic cavities might not be tuned 
correctly for the IT mode. Examination of the test hardware drawings revealed that the cavity 
configurations were different from those intended (see Figure 54). The cavity depths were too 
shallow to be optimally effective for the IT mode. Calculation using the HIFI module in 
ROCCID showed that the effects of the as tested cavity configurations on the IT chamber 
responses were indeed small (see Figure 64). 
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Effects of Cavity Configurations on Chamber Response 



Figure 64. Calculations Show Little Effect of Actual Cavity 
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IV. G, Validation Testing, (cont) 


Figure 65 is a time series showing the oscillation component of the chamber pressure 
measured in a typical unstable-combustion test. The figure shows the exponential growth of the 
amplitude typical in a linear instability. The evolutions of the amplitude and frequency of the 
chamber pressure are shown in Figure 66. The power spectral density (PSD) plot, obtained for » 

the 30 msec period during which the amplitude of the oscillations shown in the time series 
appear to be steady, is shown in Figure 67. The PSD plot in Figure 67 shows that the dominant 
frequencies were approximately 3303, 6548, and 9880 Hz. The 3303 Hz corresponds to the 
natural frequency of the first tangential (IT) acoustic mode of the chamber. The 6548 and 9880 
Hz frequencies appeared to be the harmonics of the IT resonant frequency although they may 
correspond to the natural frequencies of some higher mixed modes. Positive identification of 
higher mixed modes is difficult because differences in frequency values between one higher 
mixed mode and the next are small. 

Comparison between the measurements of chamber pressure, RP-1 and oxygen 
manifold pressures, and accelerations shows that the high-frequency pressure transducers in the 
propellant manifolds and the accelerometers mounted on the outside of the combustion chamber 
did detect instabilities and accurately provided the values of the resonant frequencies. These 
instruments were demonstrated to be well-suited in supplementing the high-frequency pressure 
transducers in the combustion chamber in combustion stability testing. 

The results shown in Figures 65 through 67 are approximately typical of all other 
unstable-combustion tests. A complete set of time series, amplitude and frequency evolution, 
and power spectral analysis of chamber pressure for all unstable tests is provided in Volume II of 
this report. 

Waterfall plots of the chamber pressures showed that instabilities began with a IT 
mode during transient start-up in tests 14, 18, 19, and 27. Approximately at the end of this start- 
up period a 2T mode popped up then disappeared immediately thereafter. The 2T frequency is 
more pronounced in the PSD plots of the chamber pressures. The waterfall and PSD plots for • 

test 14 are shown in Figures 68 and 69 as an example. The mean chamber pressure began 
decaying also at this time indicating the test was being shutdown. It is speculated that if the test 
shutdown had been delayed a little later, the instability would have switched to the 2T 
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Figure 65. Oscillation Component of Chamber Pressure Measured by PCHF1 
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Figure 66. Evolution of the Amplitude and Frequency of PCHF1 Measured 
Chamber Pressure 
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Figure 67. Power Spectral Density Analysis of PCHF1 Measured Chamber 
Pressure 
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AUTOSPECTRAL OENS1TY OF PCHF 1 (PSH 
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IV. G, Validation Testing, (corn) 


mode similar to the test 40 discussed later. It should be mentioned that all spontaneous unstable- 
combustion tests were shutdown by the combustion stability monitor (CSM) immediately after 
an instability was detected in order to preclude damage to the hardware. 

In test 17, the combustion was stable at the steady-state operating conditions. A 
large amplitude of approximately 600 psid peak-to-peak at the 1L frequency of approximately 
1300 Hz, however, was seen during early start-up transient period. Figures 70 and 71 show the 
time series and PSD plots of the chamber pressure during that time. 

In test 21, the combustion was stable during start-up transient and approximately 500 
msec of steady-state operation but became unstable during shutdown transient immediately after 
the test termination signal (FS-2). The instability began with a 2T mode then immediately 
switched to a IT mode as shown in the waterfall plot of PCHF1 (Figure 72). The PSD plot 
obtained for the period of time at the beginning of the instability clearly indicated the 2T 
frequency as shown in Figure 73. The PSD plot for a later period of time is similar to those from 
all other tests where combustion instabilities initiated during transient start-ups. It should be 
mentioned here that the injection temperature of oxygen was approximately 35 degrees higher in 
this test than those in all other tests. The warmer injection temperature of the oxygen appeared to 
have a stabilizing effect on the combustion. 

In test 29, the chamber pressure oscillated with a small amplitude (with no apparent 
growth in amplitude) during a steady-state period of approximately 0.3 sec and then suddenly the 
amplitude grew larger. 

In test 37, only PCHF1 indicated chamber pressure oscillation of small amplitude 
while all other transducers including the accelerometers were dubiously quiet. In this test, a 
large pressure spike that occurred during transient operation might have rendered a few or all of 
the high-frequency pressure transducers - and also the accelerometers although unlikely - 
inoperative or distorted. The oscillation measured by PCHF1 consisted of two major sinusoidal 
components whose peak-to-peak amplitudes and frequencies were approximately 1 15 psid at 
3008 Hz and 103 psid at 4424 Hz. The individual amplitude of each component is less than 10 
percent of the mean chamber pressure. The 10 percent value is the threshold above which engine 
operation is considered unstable as defined by CPIA publication 247. Therefore, this test 
appeared to exhibit marginal combustion stability. 
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Figure 70. Chamber Pressure Oscillation With 1 L Mode During Transient 
Start-Up of Test 17 
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Figure 71. Power Spectral Density Analysis of Chamber Pressure During 

Transient Start-Up of Test 17 Showed 1L Resonant Frequency of 
Approximately 1300 Hz 
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Figure 73. Power Spectral Analysis for the Early Part of Instability Indicated 
2T Resonance 


IV. G, Validation Testing, (cont) 


In test 39, the lowest resonant frequency was 2065 Hz which is significantly lower 
than the value of 3303 Hz commonly seen in other unstable-combustion tests. Since this fre- 
quency is not near any expected resonant mode its source is unexplained. 

Test 40 is unique in that combustion instability lasted for approximately 600 msec 
without the test being shutdown by the CSM. The CSM did not shutdown the test because the 
chamber pressure transducer that provides instability information to the CSM was inoperative. 

In fact, all chamber pressure high-frequency transducers except PCHF1 had been damaged 
during this test or during earlier tests. The long duration instability in this test provided valuable 
data that would not be known from other tests. Figure 74 is a waterfall plot of the chamber 
pressure measured by PCHF1 for the entire period from the beginning of the mean chamber 
pressure rise to the shutdown time. A higher resolution waterfall plot for only the start-up 
transient period shown in Figure 75. Figure 76 shows the mean chamber pressure versus time 
from FS-1 for the purpose of correlating the combustion instability events shown in Figures 74 
and 75 with the mean chamber pressure. Figure 74 shows that the instability began with a IT 
mode at FS 1+1.42 seconds which is during the transient start-up period (see Figure 76), and then 
switched to a 2T mode at FS 1+1. 45, a time that corresponds to the end of the transient start-up 
and the beginning of steady-state operation. The PSD analysis of PCHF1 obtained for this time 
period is shown in Figure 77. The figure clearly shows the dominant 2T frequency of 
approximately 5220 Hz. Figure 74 shows that the instability switched back to the IT mode at 
FS1+1.70 sec, which is in the middle of the steady-state operation, as noted in Figure 76. At this 
time, the operating parameters such as chamber pressure and propellant flow rates were steady. 
The cavity sound speed, however, may still have been in transient as indicated by the time series 
of the cavity gas temperatures shown in Figure 76. 

The fact that the instability mode switched from IT to 2T at the beginning of the steady-state 
operation is worthwhile to emphasize. Had the test been shutdown by the CSM before achieving 
steady-state operation like many other tests, the instability in the 2T mode would not have been 
observed. More importantly, it points out the possibility that the 2T mode might have occurred 
in several other unstable-combustion tests if they had not been shutdown so quickly by the CSM. 
This is further supported by the data from test 21 where an instability was initiated during 
steady-state operation where it began with a 2T mode. A mean chamber pressure vs time plot for 
test 1 9, which is typical of unstable -combustion tests where the test shutdowns were initiated 
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Figure 74. Evolution of the Amplitude and Frequency of Chamber Pressure in 
Test 40 Shows the Instability Began With a IT Mode, Switched to a 
2T Mode, Then Switched Back to the IT Mode 
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Figure 76. Mean Chamber Pressure and Cavity Temperature vs Time in Test 40 
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IV. G, Validation Testing, (cont) 


by the CSM, is shown in Figure 78. It can be seen from the comparison between Figures 76 and 
78 that there was little or no time for the 2T mode to develop in the unstable-combustion tests 
where the test shutdown was initiated by the CSM before steady-state operation was achieved. 
This is probably the reason for the extremely brief appearance or even non-existence of the 2T 
mode in the waterfall plots for other unstable combustion tests. 

The observed combustion instabilities have been correlated with chamber pressure, 
mixture ratio, and injection pressure drop. The stability test results are shown as a function of 
chamber pressure and mixture ratio in Figure 79. The figure shows distinct areas of stable and 
unstable operation. The unstable region is in the area of high chamber pressure and high mixture 
ratio whereas the stable region is in the area of low chamber pressure and low mixture ratio. It 
should be noted that test 37, which could not be positively characterized as previously discussed, 
was considered to be unstable in this display. Although pressure oscillation amplitudes are 
smaller in tests with monotuned and bituned cavities than those with blocked acoustic cavities, 
influences of cavity configurations could not be discerned. Again, the acoustic cavities were not 
optimally tuned because of the previously described design/drafting error. A neutral line is 
estimated on Figure 79 to define the boundary of stable and unstable regions. A precise 
determination of the stable-unstable boundary line would require more test data. The neutral line 
shows that this particular test engine may operate stably at mixture ratio less than approximately 
1.6 or at chamber pressures less than 650 psia. Operation at nominal chamber pressure of 1200 
psia and mixture ratio of 2.80 was unstable and additional damping would be required to operate 
stably. 


The stability test results were also correlated as a function of the fuel (DPF) and 
oxidizer (DPO) injector pressure drop as shown in Figure 80. This correlation also shows 
distinct stable and unstable regions. The unstable region is in the area of high DPO and high 
DPF. Again, no influences of cavity configurations are apparent from this correlation. An 
approximated stable-unstable boundary line is drawn on the figure. This line shows that this 
particular test engine may operate stably if the RP-1 injection pressure drop is less than 70 psi or 
if the oxygen injection pressure drop is less than 1 15 psi. 

H. TEST ANALYSIS AND ROCCID METHODOLOGY ASSESSMENT 

Predictions of the validation test hardware combustion stability and performance 
over a wide range of operating chamber pressure and mixture ratio were made well in advance of 
the validation testing. These predictions are described in Section IV. C. and IV.D. of this report 
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Figure 78. Mean Chamber Pressure vs Time in Test 19, Typical of Unstable- 
Combustion Tests Which Were Shutdown Early by CSM 


DPO (paid) Chamber Pressure (pai») 


open symbol : stable 

solid symbol: unstable 

number 1-20: test with bituned cavities 

number 21-28: test with blocked cavities 

number 29-40: test with a mo no tuned cavity 



Figure 79. Pc-MR Stability Map 



Figure 80. DOP-DFP Stability Map 
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IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 

and in great detail in references 14, 15, and 18. In addition, in order to provide a direct quantita- 
tive comparison of the ROCCID calculated stability and performance parameters with test mea- 
sured and derived parameters, the ROCCID methodology was applied using the point analysis 
option for each valid test. The actual hardware design parameters, including the erroneously 
designed acoustic cavities, and the measured mixture ratio, chamber pressure, propellant tem- 
peratures, and chamber throat area for each test were input to ROCCID to facilitate this direct 
comparison. 


In performing these calculations, no attempt was made to anchor any input param- 
eters using the validation test data. Anchoring parameters from Reference 13 (Table 8 of 
Reference 13) and identical to those used for the pretest predictions were used in these 
calculations. These included a factor of 1 .20 on the oxidizer dropsize, 0.08 on the atomization 
lengths for vaporization calculations, 0.9 for the predicted chamber sound speed and 0.925 on the 
overall Em predicted using Nurick's technique. In addition, the unielement Em (Em un j) was 
determined from the correlation of cold flow test data of Reference 13. 

Em un i = 1.04 - 0.185*R + 0.042*R2 - 3.37E-3*R 3 

where R is the oxidizer-to-fuel momentum ratio. The concentrated combustion plane used in 
HIFI was assumed to be at an axial position corresponding to 80% of the total energy release. 

The speed of sound in the acoustic cavity was assumed equal to that in the combustion chamber 
since the mistakenly designed cavities had large entrance openings and shallow depths. The 
validity of this assumption was not extensively investigated since the test results indicated that 
no significant damping was obtained with the cavities installed. This assumption can be viewed 
as indicative of an a-priori estimate for the purposes of the ROCCID calculations. 

A comparison of the predicted and measured performance in terms of C* and the 
thrust based energy release efficiency (ERE) is shown in Table 28. The calculated and measured 
C* combustion efficiencies are plotted against mixture ratio in Figure 8 1 . The calculated and 
measured ERE efficiencies are plotted against mixture ratio in Figure 82. It should be noted that 
all data even from tests with short duration or derived measurements were plotted in the figures. 
Both calculation and measurement showed consistencies between C* and ERE efficiencies. The 
measured performance of the engine was, however, significantly lower than predicted by the 
ROCCID methodology. 
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TABLE 28. MEASURED AND PREDICTED VALIDATION ENGINE PERFORMANCE 



Calculated 

Measured 

Test 

Pc 

Mr 

Eta-C* 

ERE 

Eta-C* 

ERE 

4 

1178 

2.59 

0.983 

0.983 

0.801 

0.842 

7 

1441 

1.45 

0.990 

0.989 

0.975 

0.972 

9 

505 

2.94 

0.988 

0.988 

0.994 

0.987 

12 

250 

1.01 

0.981 

0.980 

0.863 

0.792 

13 

1210 

1.13 

1.003 

1.004 

1.049 

1.010 

14 

1460 

2.67 

0.985 

0.985 

0.932 

0.948 

15 

969 

1.93 

0.980 

0.981 

0.968 

0.951 

17 

795 

6.74 

0.991 

0.991 

1.025 

1.015 

18 

1165 

2.35 

0.979 

0.979 

0.978 

0.997 

19 

1220 

2.06 

0.979 

0.980 

0.994 

1.000 

20 

1193 

1.50 

0.983 

0.982 

0.955 

0.948 

21 

1222 

2.67 

0.984 

0.984 

0.934 

0.928 

22 

1092 

2.75 

0.980 

0.980 

0.836 

0.924 

23 

1004 

1.25 

1.014 

1.016 

0.961 

0.937 

24 

788 

6.71 

0.990 

0.990 

0.956 

0.947 

25 

1410 

2.86 

0.987 

0.987 

0.908 

0.926 

26 

1130 

1.70 

0.983 

0.983 

0.736 

0.753 

27 

1755 

3.03 

0.990 

0.990 

0.890 

0.938 

28 

1260 

5.15 

0.994 

0.994 

0.904 

0.946 

29 

1208 

2.55 

0.982 

0.982 

0.927 

0.915 

30 

791 

3.09 

0.991 

0.990 

0.919 

0.921 

31 

953 

2.09 

0.978 

0.979 

0.843 

0.663 

36 

249 

3.31 

0.990 

0.990 

0.931 

0.899 

37 

1738 

3.05 

0.990 

0.990 

0.918 

0.927 

38 

1467 

3.07 

0.991 

0.990 

0.920 

0.963 

39 

1214 

5.60 

0.994 

0.994 

0.945 

0.916 

40 

1180 

2.23 

0.978 

0.978 

0.874 

0.886 


Note: Efficiencies greater than 1.00 are suspect and are probably the result of interpolation 

accuracies for calculated values or test measurement accuracies for measured values. 
See Section IV H 
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O Calculated ETA-C* from Table 28 
□ Measured ETA-C* from Table 28 


1.1 


See Table 25 for notes regarding data reduction. 
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Figure 81. Comparison of Unanchored ROCCID Predictions and 
Test Results - C* Efficiency 
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Isp- Based ERE 


O Calculated Igp - Based ERE from Table 28 
□ Measured l S p - Based ERE from Table 28 


See Table 25 for notes regarding data reduction. 



Mixture Ratio (0/F) 


M4.D15.16 


Figure 82 . Comparison of Unanchored ROCCID Predictions 
and Test Results - ISP-Based ERE 
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IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 


The C* efficiency from the ROCCID point analysis ranged from a low value of 
97.8% at a mixture ratio of 2.2 to a value slightly over 100% at a mixture ratio of 1.2. 
Efficiencies slightly greater than 100% as calculated by ROCCID resulted because of inflections 
in the Isp (ODE) versus mixture ratio curve as interpolated by ROCCID. No attempt was made 
to determine the vaidity of these inflection points. The ROCCID calculated ERE values were 
within 0.2% of the ROCCID calculated C* efficiencies. The measured C* efficiencies ranged 
from a low value of 87% at a mixture ratio of 2.9 to a value of 104.9% at a mixture ratio of 1.13. 
In tests with sufficient steady-state duration where accurate measurements were made, the 
measured C* and ERE efficiencies agreed within 2%. The largest difference between calculated 
and measured efficiency values is approximately 0.07 (7%). Test results showed that the 
efficiencies are more sensitive to mixture ratio than predicted. In addition, calculations showed 
that the minimum efficiency occurs at a mixture ratio of 2.2 while test results showed that the 
minimum efficiency occurs at a mixture ratio of 2.9. The disagreement suggested that the 
propellant mixing and vaporization rate might have been overpredicted. 

The vaporization efficiency at near the stoichiometric mixture ratio appears reason- 
able based on a comparison between the predicted and measured static pressure profiles in the 
combustion chamber. A typical pressure profile comparison at near nominal conditions is shown 
on Figure 83. However, at very low or high mixture ratios the ROCCID methodology overpre- 
dicts the vaporization rate, as can be seen from the comparison between the predicted and mea- 
sured chamber static pressure profiles shown on Figure 84. This overprediction is probably a 
result of the generalized length correlation used for vaporization estimates that were based on 
data and calculations assuming typical stoichiometric combustion temperatures. The ROCCID 
methodology should be modified to reflect the lower vaporization rates at off stoichiometric 
conditions either by modifying the generalized length correlation or by including a combustion 
temperature correction factor. Although Figure 84 indicates that the actual vaporization rates 
were lower than predicted in tests with low or high mixture ratios, the chamber is of sufficient 
length to produce a high overall performance efficiency. 

The disagreement between calculated and measured performance is attributed pri- 
marily to the overprediction of mixing efficiency by the ROCCID methodology for the following 
reasons: 1) The performance efficiency and therefore the vaporization efficiency is high at low 
and high mixture ratios. The efficiencies were also high at near the nominal mixture ratio 
because it can be inferred from the agreement between the measured and calculated pressure 
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Static Pressure Profile Comparison Static Pressure Profile Comparison Static Pressure Profile Comparison 

Test 29 Test 21 Test 30 

0/F=2.55 0/F=2.67 O/F=3.09 
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Figure 84. Poor Agreement Observed Between Predicted and Measured Chamber 
Static Pressure Profiles at Very Low and High Mixture Ratios 


IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 


profiles and 2) Poor mixing can result in low efficiencies when the combustor is operated near 

the nominal mixture ratio where the values of C* and ISP are maximum, but not necessarily 

result in low efficiencies when the engine is operated at very low or very high mixture ratio, 

since performance loss in fuel lean or fuel rich stream tubes is compensated for by performance y 

gains in stream tubes operating near to stoichiometric values. 

The mixing efficiency used in the performance prediction was 0.874. This value was 
obtained using the DPOPMIX correlation in ROCCID with the unielement E m of 0.77 as an 
input to the correlation and an anchoring factor of 0.93. The values of the unielement E m and an- 
choring factor were obtained from the Reference 13. A better match of the prediction to the data 
is obtained if a multielement Em value of 0.75 was used in ROCCID. 

The significantly finer pattern on the ROCCID engine compared to the Ref. 13 
engine actually decreased propellant mixing and performance, although performance was pre- 
dicted to increase with the finer pattern. Cold flow testing, shown in the photograph of Figure 57 
indicates that the propellant spray fan formation from individual elements is not complete before 
adjacent elements begin to restrict the fans. This may cause a multielement mixing efficiency 
that is actually lower than the unielement mixing efficiency. To confirm this it is recommended 
that both unielement and multielement cold flow tests be conducted to measure the cold flow 
mixing efficiencies. 


The spray fan formation observed from the ROCCID cold flow was used to estimate 
the amount of fan interaction and depict this fan overlap on the ROCCID injector pattern 
drawing. In addition, the Reference 13 injector pattern drawing was used to layout the estimated 
fan overlap from that design, using the orifice diameter ratios for the two injectors as the scaling 
parameter. The results are shown on Figure 85. The Reference 13 design shows good individual 
element fan development before intersection with adjacent fans, indicating that the presence of 
the adjacent fans do not restrict the fan formation and intra-element mixing. Of course, if the 
elements are too widely spaced, there would be no enhancement of the mixing efficiency going #, 

from a unielement to multielement configuration. The correlation currently in the DROPMIX 

* 

code indicates an increase in multielement Em for a given unielement Em as the pattern fineness 
increases. This correlation should be used with caution for very fine patterns and an expanded 
element mixing data base covering a wide range of element pattern densities and element types 
using both cold flow and hot fire data should be generated. 
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Figure 85. Comparison of ROCCID and Injector Characterization Program 
Engine Propellant Spray Fan Formation 



IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 


Nonacoustic Stability 

ROCCID calculations showed that the validation engine was chug stable at all 
operating conditions tested. This is in good agreement with test results. No chug instabilities 
were encountered in any of the tests, which covered a wide range of operating conditions. 
Calculated chug results are shown in Table 29. The calculated marginal chug PCs are plotted on 
Figure 86 for comparison with test results. Thus the ROCCID methodology is adequate in terms 
of defining injection pressure drop requirements. 

Acoustic Stability 

Acoustic stability calculated results for each test are summarized in Table 30. Test 
results are also shown in the table for comparison. It should be noted that the value of -2100 1/s 
was used for growth coefficients where the coefficients were not calculated by ROCCID because 
the mode is so stable. 

The eleven Block 1 tests with a bituned acoustic cavity covered a range of chamber 
pressures from 250 to 1441 psia and mixture ratios from 0.54 to 6.74. Seven of these tests were 
stable and four had spontaneous first tangential mode instabilities that ranged from peak-to-peak 
amplitudes of 46 to 103% of the steady-state chamber pressure. Four of these stable tests had the 
combustion process disturbed by bombs with overpressures ranging from 5 to 67% of the mean 
chamber pressure. Two stable tests were low chamber pressure tests for chug stability evaluation 
and were not bombed and one test had a bomb intended but did not produce a measurable over- 
pressure. The ROCCID methodology correctly predicted all of the four unstable tests but only 
two of the seven stable tests. 

The eight Block 2 tests without an acoustic cavity had chamber pressures from 788 
to 1706 psia and mixture ratios from 1.25 to 6.71. Only two of these tests were stable, the lowest 
and highest mixture ratio tests. These tests had bomb overpressures of 5 and 1 1% of the mean 
chamber pressure. Of the six unstable tests all were in the first tangential mode and all were cor- 
rectly predicted to be unstable by the ROCCID methodology. The methodology correctly pre- 
dicted, however, only one of the two stable tests. 

The eight Block 3 tests with a monotuned cavity had chamber pressures from 249 to 
1735 psia and mixture ratios from 2.09 to 5.60. Two of these tests were stable and six were 
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Table 29. ROCCID Calculated Chug Results 
ROCCID Version 30-Jul-91 
Date: 23-Oct-91 


| Calculated Results | 


+- 








1 

Test 

Pc 

MR 

i 

Marginal 

F 

i 

1 

1 


(psia) 


i 

i 

Chug Pc 
(psia) 

(HZ) 

i 

i 

1 




i 



i 

1 

i 

4 

1178 

2.59 

i 

i 

302 

705 

i 

i 

1 

7 

1441 

1.45 

i 

258 

608 

i 

1 

9 

505 

2.94 

i 

257 

645 

i 

1 

12 

250 

1.01 

i 

83 

744 

i 

1 

13 

1210 

1.13 

i 

117 

763 

i 

1 

14 

1460 

2.67 

i 

288 

694 

i 

1 

15 

969 

1.93 

i 

280 

662 

i 

1 

17 

795 

6.74 

i 

76 

658 

i 

1 

18 

1165 

2.35 

i 

291 

696 

i 

1 

19 

1220 

2.06 

i 

303 

704 

i 

1 

20 

1193 

1.50 

i 

171 

744 

i 

1 

21 

1222 

2.67 

i 

274 

670 

i 

1 

22 

1092 

2.52 

i 

303 

710 

i 

1 

23 

1004 

1.25 

i 

142 

759 

i 

1 

24 

788 

6.71 

i 

82 

681 

i 

1 

25 

1410 

2.86 

i 

277 

689 

i 

1 

26 

1130 

1.70 

i 

269 

633 

i 

1 

27 

1755 

3.03 

i 

281 

684 

i 

1 

28 

1260 

5.15 

i 

118 

683 

i 

1 

29 

1208 

2.55 

i 

282 

676 

i 

1 

30 

791 

3.09 

i 

258 

654 

i 

1 

31 

953 

2.09 

i 

281 

670 

i 

1 

36 

249 

3.31 

i 

224 

600 

i 

1 

37 

1738 

3.05 

i 

257 

657 

i 

1 

38 

1467 

3.07 

i 

260 

671 

i 

1 

39 

1214 

5.60 

i 

85 

651 

i 

1 

+-• 

40 

1180 

2.23 

i 

274 

670 
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Comparison of Unanchored ROCCID Predictions and Test Results 

Marginal Chug Pc 



Figure 86. Comparison of Unanchored ROCCID Predictions and Test Results 
Marginal Chug Pc 
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TABLE 30. 

WCCID Version 30-JU.-91 COMPARISON BETWEEN CALCULATED RESULTS AND TEST DATA 

Oste: O0*OCT*91 
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[1) Chamber pressure below which chug instability is predicted for given test mixture ratio. 

[2) AL is the growth coefficient calculated by ROCCID as illustrated in Figure 18. 

[3) Calculated results using ROCCID as described in Section IV H. 


IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 


unstable. Similar to the Block 2 comparison results, the ROCCID methodology predicted cor- 
rectly all of the six unstable tests but only one of the two stable tests. 

The combustion stability test result from all three test blocks are shown in Figure 87 
as a function of measured chamber pressure and mixture ratio. An approximate correlation for 
the stable/unstable transition as a function of chamber pressure and mixture ratio is also provided 
and compared to the ROCCID prediction. From these results, we were not able to discern any 
change in this transition boundary with the use of the improperly tuned acoustic cavities. Some 
additional damping was evident when these acoustic cavities were used, however, based on the 
observed reduction in the amplitude of the instability and reduced experimental growth coeffi- 
cient. Similar trends were predicted by ROCCID as noted in Table 30. 

A comparison of the experimentally observed stability transition boundary (from 
Figure 87) to that predicted by ROCCID (IT lower boundary from Figure 35) is shown in Figure 
88. From this comparison, the zone of unstable operation expected based on the ROCCID results 
was somewhat larger than the experimentally determined zone. Interestingly, however, a close 
correlation between the ROCCID prediction and the observed test results was achieved near the 
design mixture ratio of 2.8. This may indicate that the stability analysis models provide better 
predictions at mixture ratios near "nominal" values where a majority of the historical database 
was generated and upon which the models were developed and verified. 

Although all unstable-combustion tests were correctly predicted to be unstable, the 
instability mode and the growth coefficients were not correctly predicted in several of the tests. 
The capability to predict precisely the instability modes and accurately the growth coefficients 
require further improvements of the component models used within ROCCID. As previously 
discussed in Section IV.G, early shutdown of the test by CSM after a IT instability was detected 
during transient start-up provided little or no time for steady-state operation at which the insta- 
bility may switch to a 2T mode similar to what was seen in test 40. This might explain the dis- 
crepancies between the calculated instability mode and test results in a few of the tests. 

The discrepancies between calculated and test results are believed to result also from 
the use of the empirical combustion response and the chamber response model. The weakness of 
the models is discussed in the following paragraph. 
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Figure 87. ROCCID Validation Test Results Indicate Stable 
and Unstable Operating Regions 


170 


150 


1 



08 S/LUq| ‘ 81 RIM 0 IJ l~dU 


171 


Than was Observed Experimentally 


IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 


The combustion response used in the predictions follow Crocco and Cheng theory, 
which uses to the pressure interaction index, n and a sensitive time lag, t. The value for n was 
. obtained using Smith and Reardon's empirical rule - by equating it to the time it takes for 20% of 
the controlling propellant to evaporate. The controlling propellant is defined as the propellant 
that is less volatile than the other. For the LOX/RP-1 propellant combination, RP-1 is the con- 
trolling propellant. The vaporization time, hence the sensitive lag, is calculated using Priem's 
generalized length model in ROCCID with the propellant dropsize calculated using Aerojet 
atomization model. Thus, the accuracy in the prediction of the combustion response depends on 
the accuracies of the atomization model and the vaporization model. 

Interestingly, as noted earlier in this section, the propellant vaporization rate appears 
to be over predicted at the very low and high mixture ratios where the combustion temperature is 
significantly lowered. The over prediction in vaporization for the extreme mixture ratios would 
lead to smaller calculated values of sensitive time lags and thus a larger calculated zone of insta- 
bility with respect to mixture ratio. 

In spite of the relatively good correlation between the ROCCID predictions and the 
validation test results, determination of the combustion response remains a highly empirical and 
analyst dependent process. The adequacy of a pirori specification of the combustion response 
can be a function of the element type and design, the propellants and the available experimental 
data base. While some analytical approaches, such as CRP and the NASA Lewis Injection 
Coupling Model have been developed, their application has been relatively limited and generally 
involved in correlation of tested designs. Clearly, further advancement in the modeling and 
characterization of the combustion response is essential for improvement in the overall predictive 
capability of liquid propellant combustion stability. 

However, since the whole sensitive time lag methodology, upon which ROCCID is 
based, is built upon empirical correlations and analytical approaches that have been verified, 
refined and/or anchored with empirical data generated over the past thirty to thirty-five years, it 
is not obvious that additional improvements in any specific modeling area would result in an 
overall improvement in the combustion stability characterization. For example, the chamber 
response modeling does not account for the acoustic energy damping provided by the chamber 
wall and the propellant spray. Improvements in this area may be able to be incorporated into the 


RPT/F0213.105-IV 


172 


9/30/92 


IV. H, Test Analysis and ROCCID Methodology Assessment, (cont) 

chamber response model to make it more physically representative. Yet since, the combustion 
response, n, has been "backed out" from experimental results using chamber response models 
without the wall and spray damping considerations, its use with the improved chamber response 
models may actually result in a degradation in the predictive capability of the methodology. 

By combining the best performance and stability models into one program and 
giving them a standard base for comparison, ROCCID has made it possible to rigorously 
evaluate the models incorporated into the program. Until better models can be developed to 
accurately and consistently predict the critical parameters that affect engine performance and 
stability, predictions by the models in ROCCID will have a large error band. Improved 
diagnostic equipment will permit the acquisition of better data to improve and validate the 
models. More mechanistic models can be incorporated into ROCCID, which require fewer 
assumptions in their calculations. CFD generated empirical models for portions of the 
combustor can also be incorporated. 
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